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PREFACE. 



In March, 1875, the usual biennial notice was issued, giving the 
subjects for the Adams Prize to be adjudged in 1877. The 
following is the chief portion of the notice : 

The University having accepted a Fund raised by several members 
of St John's College for the purpose of founding a Prize to be called the 
Adams Frize, for the best essay on some subject of Piu*e Mathematics, 
Astronomy or other branch of Natural Philosophy, the Prize to be 
given once in two years, and to be open to the competition of all 
persons who have at any time been admitted to a degree in this 
University — 

The Examiners give notice that the following is the subject of the 
Prize to be adjudged in 1877 : The Criterion qfJh^namical StahiMy. 

To illustrate the meaning of the question imagine a particle to slide 
down inside a smooth inclined cylinder along the lowest generating line, 
or to slide down outside along the highest generating line. In the 
former case a slight derangement of the motion would merely cause 
the particle to oscillate about the generating line, while in the latter 
case the particle would depart from the generating line altogether. 
The motion in the former case would be, in the sense of the question, 
stable, in the latter unstable. 

The criterion of the stability of the equilibrium of a system is, 
that its potential energy should be a minimum; what is desired is, a 

R. A. h 
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corresponding condition enabling us to decide when a dynamically pos- 
sible motion of a system is such, that if slightly deranged the motion 
shall continue to be only slightly departed from. 

The essays must be sent in to the Vice-Chancellor on or before the 
16th December 1876, &c., &c. 

S. G. PHEAR, Vice-chancellor. 

J. CHALLIS. 

G. G. STOKES. 

J. CLERK MAXWELL. 



The pressure of other engagements for some time prevented 
me from giving my attention to the subject. This essay was 
therefore almost entirely composed during the year 1876. It is 
now printed as it was sent in to the Examiners, the changes being 
merely verbal. Some few additions have been made where ex- 
planation appeared to be necessary, but all these have been 
marked by square brackets, so that they can be at once dis- 
tinguished from the original parts of the essay. 

In order to shorten the essay as much as possible many 
merely algebraic processes have been omitted and the results only 
are stated. It is hoped that this will add clearness as well as 
brevity to the reasoning, as the attention of the reader will not 
be called from the argument to follow a manipulation of symbols 
which may not present any novelty. 

The line of argument taken may be indicated in a general way 
as follows. Chapter i. begins with some definitions of the terms 
stable and steady motions. It is then pointed out that whether 
the forces which act on the system admit of a force-function or not, 
the stability of the motion, if steady, is indicated by the nature 
of the roots of a certain determinantal equation. The boundary 
between stability and instability being generally indicated by the 
presence of equal roots, a criterion is investigated to determine 
beforehand whether equal roots do or do not imply instability. 
This case being disposed of, the consideration of the determinantal 
equation is resumed. Two general methods are given by which, 
without solving the equation, it may be ascertained whether the 
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character of the roots Imply stability or instability. These occupy 
Chapters ii. and ill. In the first method a derived equation is 
made use of, and it is shown that a simple inspection of the signs 
of the coeflScients of the several powers in these two equations will 
decide the question of stability. In the second method a certain , 
easy process is found which if performed on the determinantal 
equation will lead to the criteria of stability. At the end of the 
third Chapter a geometrical interpretation is given to the argu- 
ment. 

In the fourth Chapter the forces which act on the system 
are supposed to have a force-function. The determinantal equa- 
tion is then much simplified. Several points are considered in 
this Chapter which are necessary to the argument, such as the 
proper method of choosing the steady co-ordinates (if there be 
any), the distinction between harmonic oscillation about steady 
motion and that about equilibrium, and the changes which must 
be made in the determinantal equation when the equations of 
Lagrange become inapplicable. A method of modifying the 
Lagrangian function is also given by which, in certain cases, 
the fundamental determinant may be reduced to one of fewer 
rows and columns. 

In the fifth Chapter a series of subsidiary determinants is 
formed, and it is shown that at least as many of the conditions 
of stability are satisfied as there are variations of signs lost in the 
series in passing from one given state to another. It is also shown 
that this is equivalent to a maximum condition of the Lagrangian 
function. 

In the sixth Chapter the energy test of stability is considered. 
It is also shown that, when the motion is steady, this reduces 
to the same criterion as that indicated in Chap. V. 

In the seventh Chapter the question considered is whether 
the stability of a state of motion can really be determined by an 
examination of the terms of the first order only. In some cases 
these are certainly suflScient, and an attempt is made to discrimi- 
nate between these cases and those in which the terms of the 
higher order ultimately alter the character of the motion. 

• If the Hamiltonian characteristic and principal functions be 
given, the conditions of stability as regards space only, or both 
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space and time may be deduced. But if these be not known as 
expressed in the Hamiltonian form, we may yet sometimes dis- 
tinguish between stability and instability if we can determine 
whether a certain integral ceases to be a minimum at some instant 
of the motion. This is the subject of the eighth chapter. 

As part of the third edition of my treatise on the Dynamics of 
Rigid Bodies was written at the same time as this essay, there are 
necessarily points of contact between the two works. Thus the 
subjects of the first part of the seventh chapter and of a portion 
of the sixth will be found discussed in the treatise on Dynamics. 
But as the objects of the two books are not the same, it will be 
found that in all these cases there are considerable dififerences in 
the modes of demonstration. 

EDWARD J. ROUTH. 



Petbbhouse, 

August 14, 1877. 
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CHAPTEE I. 



DeJiniUona of the terms small quantity, stable motion, steady motion. 
Arts. 1, 2. 

A system of bodies in steady motion is stable if the roots of a certain 
determinomtdl equaiion a/re such that their reed pa/rts a/re aU negative. 
Art 3. 

affect of equal roots, and a test to determine whether equal roots do or do 
not introduce terms which contain the tims as a factor. Arts, 4 — 7. 

Object of Chapters II, and III, Art, 8. 

1. Let us suppose a dynamical system to be set in motion 
under any forces and to move in some known manner. If any 
small disturbance be given to the system, it may deviate only 
slightly from its known motion, or it may diverge further and 
further from it. Let 0, (f>, &c. be the independent variables or 
co-ordinates which determine the position of the system, and let 
the known motion be given by ^ = ^^j ^ = ^o» &^- where ^^, <f>Q, &c. 
are known functions of the time t. To discover the disturbance of 
the system we put = 0q + x, <f>=(f>Q-^y, &c. These quantities 
X, y, &c. are in the first instance very small because the disturb- 
ance is small. The quantities x, y, z, &c. are said to be small 
when it is possible to choose some quantity numerically greater 
than all of them, which is such that its square can be neglected. 
This quantity may be called the standard of reference for small 
quantities. 

If, after the disturbance, the co-ordinates x, y, «, &c. remain 
always small, the undisturbed motion is said to be stable; if, on the 
other hand, any one of the co-ordinates become large, the motion is 
called unstable. 

It is clear that the same motion may be stable" for one kind of 
disturbance and unstable for another. But it is usual to suppose 
the disturbance general, so that if the motion can be made un- 
stable by any kind of disturbance (provided it be small) it is said 
to be unstable. On the other hand, it will be called stable only 
wlien it is stable for all kinds of small disturbances. 

R. A. 1 
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2 DEFINITIONS. [CHAP. 

2. To determiDe whether x, y, z, &c. remain small, we must 
substitute for 0, <f>, &c. in the equations of motion their values 
^o + ^» 0o + y> <^c« Assuming that x, y, &c. remain small, we may 
neglect their squares, and thus the resulting equations will be 

linear in x, y, z, &c. The coefficients of a?, -^ , ^ > y> ^ > j^> &c- 

in these equations may be either constants or functions of the 
time. In the former case the undisturbed motion is said to be 
steady for these co-ordinates, in the latter tmsteady. In the case 
of a steady motion x, y, z, &c. are all functions of the time 
which has elapsed since the disturbance and of certain constants 
of integration which are determined by the initial values of 

€uB du 

fl?, ^ , y, ^ , &c. We may therefore define a steady motion to be 

such that the same change of motion follows from the same initial 
disturbance at whatever instant the disturbance is communicated 
to the system. 

If all the coefficients in the equations to .find x, y, z are con- 
stant, they may be made to contain ^ by a change of co-ordinates. 
Thus we may write for a,, y, «, &c. 

«?=af + )8i7 + ... 

y=a'f + )8'i7+... 

z = &c. 

where a, /8, &c. are any functions of t we please. Conversely, when 
the coefficients are functions of t, we may sometimes make the 
coefficients constant by a proper change of co-ordinates. But this 
cannot always be done. If there are n co-ordinates, we have w* 
arbitrary functions a, )8, &c. at our disposal. In each of the n 
linear equations of motion we may have three terms for each 
co-ordinate, and thus we have (3 w — 1) w coefficients to make con- 
stants. We have therefore in general too many equations to 
satisfy. The proper method of choosing the co-ordinates of refer- 
encQ will be considered in a future chapter. 

3. Let us suppose a dynamical system to be making small 
oscillations under the action of any forces which may, or may not, 
possess a force function and to be subject to any resistances which 
vary as the velocities of the parts resisted. The general equations 
of motion will then be of the form 



(^''S+^''l+<)«'+(^»'S+^''l-^^»>+^°-=^ 



&c. = 
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I.] THE DETERMINANTAL EQUATION. 3 

To solve these equations we write 

x^Me"^, y^M€^, &c. 

Substituting we obtain a determinantal equation to find m. 
If we put 

A! = A^ffif + Aim + A^, &c. 

this equation may be written in the simple form 



A, B, C, 

A\ ff, C\ 



= 0. 



We may also write the equation in the form/(w) = 0. 

The coefficients M, M\ &c. are not independent, but if we 
represent the minors of -4, 5, (7, &c. by a, J, c, &c. we may easily 
show that 



We also have 



M M M" , 
— = ^- = — = &c. 
a 






It may be shown by properties of determinants that these 
equations all give the same ratios. If A, be the second minor 
obtained from the determinant /(m) by omitting the first and 
second rows and columns^ we know that 

A,/ (m) = a5' - dh. 
Hence if /(w) = we have - = t? . 

Qi 

In the same way we may show that -> = - , and so on. This 

Cb C 

property of Determinants is given in Dr Salmon's Higher Algebra, 
Lesson iv. Ex. 1. 

The general solution of the equation may therefore be written 
in the form 

z = &c. i 

where Zj, i, ... are arbitrary constants, a^, a,, &c. the values of 
the minor a when m^, m, . . . are substituted for m ; 6^, 6^ . . . the 

1—2 



4 NATURE OF THE STABILITY [CHAP. 

values of the minor b when similar substitutions are made, and 
so on. 

4. We see that the whole character of the motion will de- 
pend on the signs of the quantities m^, m,... If any one be real 
and positive, x, y, &c. or some of them will ultimately become 
large, and the steady motion about which the system is oscillating 
wiU be unstable. If all the roots are real, negative or zero and 
unequal, the motion will be stable. 

If two of the roots be imaginary we have a pair of imaginary 
exponentials. If these imaginary roots be a + $aJ — 1, the terms 
can be rationalized into 

e^{N^ cos I3t + N^ sin ^t). 

The motion will be stable if a be negative or zero, and unstable 
if a be positive. 

If two roots be equal, the form of the solution is changed. Let 
m^ = m^ + h where h will be ultimately zero, we then have 

If we now make L. and i, infinite in the usual manner, we 
find 






&c. = &c., 
where Jf^, M^ are two arbitrary constants which replace L^y L^. 

In the same way if three roots are equal we have 

This rule will be found convenient in practice to supply the 
defect in the number of arbitrary constants produced by equal 
roots. At present we are only concerned with their eflfect on the 
stability of the system. The terms which contain ^ as a factor 
will at first increase with ^, but if m be negative, the term t^e^ 



can never be numerically greater than — . If w be very small 



n 



I.] IS INDICATED BY THE ROOTS. 5 

the initial increase of the terms may make the values of x and y 
become large, and the motion cannot be regarded as a small 
oscillation. But if the system be not so much disturbed that 



em 



is large, the terms will ultimately disappear and the 

motion may be regarded as stable. If, however, the real parts of 
the equal roots are positive or zero, the terms will become large 
and the motion will be unstable. 

5. In some cases, however, the relations which exist between 
the coefficients are such that the terms which contain ^ as a factor 
are all zero. It is of some importance to discriminate these cases, 
for the stability of the system is then imaffected by the presence 
of equal roots. 

Let us suppose first that the determinantal equation has two 
roots only equal to m^, and let the terms depending on these be 



c. = &c. 



II. 



&c. = &c. 

Substituting in the equations of Art. (3) we have, following 
the same notation as before, 

an^+bn; + gn;' + ... = o^ 
^'iv;4-5'Jv;+c7xv;'+... = oi i. 

&c. = J 

&c. = &c. 

To avoid entering more minutely than is necessary into the 
properties of linear equations, we shall assume that these equations 
for the given value of m lead to but one solution with two of the 
-ATs arbitrary, unless the determinantal equation has more than 
two roots equal to m^. If in this unique solution the N^& are dll 
zero we must have 

^'i\r-h^iV;H-... = o[ III. 

&c. =0 



6 EFFECT OF EQUAL BOOTS [CHAP. 

Since two of the constants N.^ N.\ &c. are to be arbitrary, let 

N N' 
them be N.. N' then since — * = -^i- we must have the minors a 
* * a 

N N'' 
and b each equal to zero. Also since — ^ = — ^, we shall have N^' 

Ott c 

infinite unless c = 0. In the same way we may prove that all the 
other first minors are zero. 

And if the first minors are zero, we may show that two of the 
equations may be deduced from the others. Let the symbol 

IAfjy] represent the second minor, with the usual sign, formed 
y omitting the rows and columns in which AB A!S occur. 

Then since the minors a, &, c, &c. are zero, we have 



^UJ"^^''U'J?''J+---=^ 



JV. 



&c. = 0. 

Omitting the first line of III. let us multiply the others by 
.,x^ .„-pit &c., respectively. Adding the results, we have 

an identity. Hence the second equation may be deduced from 
the others which follow it. In the same way, the first equation 
may be deduced from the others. 

Rejecting the first two equations, let us transpose the arbitrary 
constants N. and Nl to the right-hand sides of the remaining 
equations. If there are to be only two arbitrary constants, these 
remaining equations must be independent; solving, we have 

with similar equations fbr the others. Hence the constants 
N^, Nl, &c., are connected by equations of the form 

^.[IS]-^/[i?]-^"[il]-<'. 

SO that when any two are chosen as the arbitrary ones, the others 
may be deduced from them. 

V 

If the determinantal equation has three roots equal to m^, and 
if the terms which contain ^ as a factor are all zero, the equa- 
tions III. must admit of a solution with three of the constants 
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N^, N^', &c. arbitrary. If these be N^^ i^T/, i^T/', we see that the 
second minors zvpr J'r" A'Ti'l ^^®* "® ^®^^* 
Next since 

we must have n/rj; ai^d j^tv both zero or ^/" infinite. 

Hence all the second minors are zero. And if the second 
minors are all zero, we have 



I// 



ABO 

A'ffC 

_A"B'G"\ 



+ A 



III 



' ABG ; 



+ ... = 0, 



and by similar reasoning three equations may be deduced from 
the remaining ones. We have then 



i^. 



" BCD ' 


-K 


' AGD ' 


+N^' 


■ ABD ' 


-Nl" 


' ABG ' 


EG'V 




A'C'jy 




A'BU 




ABG' 


B'C'U'. 




lA" c"iy'\ 




A'B'B', 




A"B' O', 



=0, 



with similar equations. 

Since f[rn) is the determinant formed by eliminating N^, N^\ 
&o, from III. we have 

dfjm) ^df{m)dA ^ df(m)dB ^ ^^ ^ df(m)dA' ^^ 
dm dA dm dB dm ' dA! dm 

dA.dB..^^ ,dA'. 



dm dm 



dm 



This vanishes when a, b, &c., a', &c. are all zero. If therefore the 
first minors of /(m) all vanish when m = m^, the equation /(m) == 

has two roots equal to m^. In the same way -v— vanishes if all its 

first minors are zero. But 

cT/ (m) _ d^A dadA o 
dm dm dm dm 

dn 

vanishes if a, ^j— , &c. are all zero. If therefore the first and 

dm 

second minors of /(m) all vanish when m^m^^ the equation 
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f[m) =0 has three roots equal to m^. It is evident the proposi- 
tion may be extended to any nuiAber of roots. 

The test that, when equal roots occur in the determinantal equa- 
tion^ the terms in tlie values of x, y, &c. which contain i as a factor 
sfumld be absent may be stated thus. If tliere are two equal roots 
all the first minors must vanish. If three equal roots, all the first 
and second minors must vanish, and so on. In these cases the 
equal roots introduce merely a corresponding indeterminateness 
into the coefficients. 

When there are more equal roots than there are rows in the 
determinantal equation, it is easy to see that there must be some 
terms in the integrals which contain ^ as a factor. 

[The following simple example will illustrate the application 
of this test. 

A particle is in equilibrium at the origin of co-ordinates under 
the action of forces whose force function U is given by 

U^^Ax^'+iBy^ + ^C^ + Dyz + JEzx + Fxy. 

If the level surfaces are ellipsoids and the force acts inwards, 
it is* clear that the equilibrium of the particle must always be 
stable. If then any equal roots occur in the determinantal 
equation, the test should show that the terms which contain t as 
a factor are absent 

If T be the semi vis viva of the particle and if its mass be 
taken as unity, we have 

Omitting accents and forming the discriminant of — wTh- U 
we have the following determinantal equation : 

= 0. 



A-m' 


F 


E 


F 


B-m' 


D 


E 


D 


C-m* 



This is the "discriminating cubic" which determines the axes 
of the quadric ?7= c, where c is a constant. The conditions that 
two of its roots should be equal, i.e. that the quadric should be 
a spheroid, are well known to be 

. EF ^ FD ^ DE 

where m^ is equal to either root. These are just the conditions 
obtained by equating any first minor of the determinsuit to zero. 
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The conditions that three of the roots of the cubic should be 
equals i.e. that the quadric should be a sphere, are 

These are the conditions that every second minor should vanish. 

In this example we have taken the case of a single particle. 
Similar remarks however apply when any system !i Indies is 
disturbed from a state of stable equilibrium. The oscillations 
may be found by the method of Lagrange. The final determi- 
nantal equation may be conveniently formed by equating to zero 
the discriminant of '-m*T+ Z7, where Tis the semi vis viva with 
the accents denoting diflferentiations with regard to the time 
omitted, and U is the force function. It is a known theorem that 
the existence of finite equal roots does not affect the stability of 
the equilibrium. Hence the conditions for equal roots must be 
such as to make all the minors equal to zero. Conversely, this 
theorem will often conveniently give the conditions that Lagrange's 
determinant has equal roots.] 

6. That there should be a difference in the modes in which 
equal roots affect the motion is no more than we should expect 
a priori. Suppose the coeflScients of the equation /(m) =0 to be 
functions of some quantity n, and that as n passes through the 
value Uq, two roots become equal to each other. Let the quadratic 
factor containing these roots be m' + 2am + )8, and let us consider 
only the case in which a and fi are real. We have a" — )8=0 
when n = w^. If a^ — )8 change sign as n passes through the value 
Wjp the roots will change trom a trigonometrical to a purely ex- 
ponential form, which would indicate a change from oscillatory to 
non-oscillatory motion. The passage from one kind of motion to 
the other may be effected through a motion represented by expres- 
sions having the time as a factor. But if a" — )8 does not change 
sign, for example, if it be a perfect square for all values of n, there 
will be no change from one kind of motion to the other, and in this 
case we should expect that the motion when the roots are equal 
will be represented by terms of the same character as before. 
Briefly, we may expect equal roots to introduce terms with ^ as a 
factor at the boundary between stability and instability ; and to 
introduce merely an indeterminateness into the coeflScients when 
the motion is stable on both sides. 

It is easy to show that in the first of these two cases the 
minors could not contain either of the factors of m" -h 2am + jS, 
For since of — jS changes sign, these factors are in one case ima- 
ginary ; and therefore if one factor occur in any minor the other 
must also be present. The minors would not only vanish, but 
must have equal roots also. But as in Art. (3), \f{m) =^ab* —ab. 
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Hence if all the first minors have equal roots it is clear that either 
/(m) has more than two equal roots, or all the second minors must 
vanish. The latter is impossible unless f{m) has more than two 
equal roots. 

These general considerations are not meant to replace the 
proofs given in the last article, but merely to explain how a 
difference in the effects of the equal roots might arise. 

7. Summing up what precedes, we see that if a dynamical 
system have n co-ordinates its stability depends on the nature of 
the roots of a certain equation of the 2n\h degree. 

If the roots of this equation are all unequal, the motion will 
be stable if the real roots and the real parts of the imaginary roots 
are all negative or zero, and unstable if any one is positive. If 
several roots are equal the motion will be stable if the real parts 
of those roots are negative and not very small, and unstable if 
they are negative and small, zero, or any positive quantity. But 
if, as often happens in dynamical problems, the terms which con- 
tain ^ as a factor are absent from the solution, the condition of 
stability is that the real roots and the real parts of the imaginary 
roots of the subsidiary equation should be negative or zero. 

8. When the equation /(Z)) = is of low dimensions we may 
solve it or otherwise determine the nature of its roots; the 
stability or instability of the system will then become known. 
But if the degree of the equation be considerable this is not a very 
easy problem. We shall devote the two next Chapters to the 
consideration of two methods by either of which, without solving 
the equation, we can determine the conditions that the real roots 
and the real parts of the imaginary roots should be all negative. 
The determination of these conditions has, it appears, never before 
been accomplished.* The consideration of the equations of motion 
will then be resumed, and the form of the determinantal equation 
f{D) = when the forces admit of a force function will be more 
particularly investigated. 

* [These conditions for the cases of a biqnadratio and a qnintic had been found 
by the Author in 1873, and read before the London Mathematical Society in June, 
187i. See also the third edition of the Author's Rigid Dynamics, Art. 436.] 



CHAPTER n. 



Statement of the theory hy which the necessary and sufficient tests of 
stability are/oimd. Objections to this theory. Arts, 1 — 3. 

These tests shown to be integrtdJmuMons of the coefficients. Art. 4. 

Method of finding these tests when the coefficients of the eqtmtian are 
mimertcal, or when several terms are absent. Arts, 5, 6. 

All th^e tests shovm to be derimUe from one called the ftmdamental 
term. Arts. 7, 8. 

The fundamental term found as an elimvrumt. Art, 9. 

A method of finding the fundamental term by derivation from the fimda- 
mental term of one degree lou)er. Art. 10, 

Another and better method of doing the same by means of a differential 
equaiwn. Arts, 11 — 13. 

1. The object of this Chapter has been explained at the end 
of Chapter I. Briefly, the criterion that the motion of a system 
of bodies should be stable is that the roots of a certain equation 
should have all their real parts negative. We propose to investi- 
gate these conditions. 

Let the equation to be considered be 

f{x) ^'p^^'Vp^oT^ + ... +jPn-i« + P»= 0. 
Let the real roots be a^> a, ... and the imaginary roots be 

a,±)8,V-l, a,±)8,\/^=l, &c. 

Then 

/(^)=i'o («-«i) (^-«i) ••• (a?"- 2a,aj + a,«+i8i'), &c. 

If then a^, a^, &c. a^, ti^^ &c. are all negative, every term in 
each factor, and therefore in the product, must be positive. 

It is therefore necessary that every term in the equation 
/(a?) = should have the same sign. It will be convenient to 
suppose this sign to be positive. 
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It is also clear on the same suppositions that none of the 
coeflScientspo»l>i' •••P» can be zero, except when the roots of the 
equation are all of the form ±)8\/— 1, or when some of the roots 
ai*e zero. 

2. Let us now form the equation whose roots are the sums of 
the roots oif{x) taken two and two. Let this be 

F{x)^P,ar-\^P^x'^'+ ... + P^^x +P.„ = 0, 

n— 1 . . . ' 

where m^n ^ - . The real roots of this equation will be 

a^ + a^, a^ + ctg, &c. 2x^, 2afj, &c. and the imaginary roots will be 

a^ + flfj + )8j V — 1, &c. It is clear from the same reasoning as 
before that if a^,a^y&G. a,f ffj, &c. are all negative, the coeflScients 
P^ Pj, &c. must all have the same sign. 

Conversely, if p^, p^ . . . have all the same sigD, the equation 

f{x) can have no real positive root, and if P,,, P^ ... P^ have all 

the same sign the equation -^(0?) can have no positive root, and 

therefore f{x) can have no imaginary root with its real part 

positive. 

3. Ovr first test of the stability of a dynamical system is that 
all the coefficients of the dyna/mical eqttation f (D) = and all the 
coejfidents of its derived equation F (D) = should have the same 
sign. 

It should be noticed that though these conditions are all 
necessary and sufficient, they are not all independent. We obtain 
too many conditions. In many cases, however, we can at once 
reduce them to the proper number of independent conditions, and 
when this is difficult we can have recourse to the second method, 
to be given in the next Chapter, which is free from this objection. 

In order to apply this method with success, it is necessary to 
have some convenient methods of calculating the coefficients 
P P P 

4 The first method which suggests itself is one similar to that 
usually given to determine the coefficients of the equation whose 
roots are the squares of the diflferences of the roots of any given 
equation. 

If £»j, £»j... be the sums of the first, second powers, &c. of the 
roots of the equation fx = 0, we have by Newton s theorem 

where p^ has been put equal to unity. If 2^ S,... be the sums of 
the powers of the equation F{x) = 0, we have in the same way 

2 4-P2 4-P2 + =0* 
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we may also prove 

2,= (n-l)af„ 

2, = (n-2)5, + S;, 

2.= (»-4)fi'.+ 3^,fl;, 

1, = in-8) 8, + 4,8,8, + Z8*, 

2, = (n - 16) 8, + 58,8, + 108,8, i 

and the general relation can be found without difSculty. 

In this way we find 
P. = (n-l)i,,. 

„_ (r.-l)(n-2)(n-3)(n-4) , («-2)'(n-3) 
*~ 1.2.3.4 P^"^ 1.2 ^»^» 

+ (n-3)>,;,, + ("-f/^"^) j>/+(n-8)p.. 

But the process becomes longer and longer at every stage. 
We shall therefore proceed to point out some other methods of 
obtaining the coefficients. 

This method of proceeding has indeed been stated only because 
it proves in a convenient way that when Pq= 1, all the coefficients 
P^, P^, Pg... of the derived equation are integral rational func- 
tions of the coefficients^^, p^... p^, 

5. The equation f (x) = beinff given, to calculate the coefficients 
o/F(x)=0. 

Put x = y ±z and equate separately to zero the sums of the 
even and odd powers of z, we have 



^ . «♦ 



/(y)+/"(y)^+r(y)g+-=o 



z* 



/'(y)^+/"'(y)^+- =0 

Rejecting the root ^ = 0, let us eliminate z. Then the roots 
of the resulting equation in y are the arithmetic means of the 
roots of/(x) = 0. 
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If, on the other hand, we eliminate y we have an equation of 
an even degree to find z. This, putting 4«* = f, is the equation 
whose roots are the squares of the diflferences of the roots of the 
given equation. 

It may be thought that this elimination may prove tedious, 
but it will be presently shown that only the first and last terms of 
the result are really wanted. All the others may be omitted in 
the process of elimination, and thus the labour will be greatly 
lessened. The method is however most useful when the given 
equation has several of its terms absent. 

6. Example, To determine the condition that the roots of the 
hiqwadraJtio 

X* +px' + qx' + rx + s = 

should indicate a stable motion. 

Applying the rule we have 
F(x) =a?' + 3pajP+ (3/ 4-2j) a?* + (4pg +/) a? 

+ (2p^q + pr + ^'-4is)a?+ {pg^ +p^r -- itps) x+pqr --r^ - p^s. 

The first four coefficients contain only positive terms, and need 
not be considered. If the last three coefficients be called P^, P^, 
P^, we have 

/P,-4P,= (|)gr-2r)» + 2i>V + pV, 
i>P,-4P, = (;)gr-2r)' + pV. 

If then Pg is positive, all the other coefficients are positive. 

The necessary and sufficient conditions of stability are there- 
fore that p, q, r, s should be finite and positive, and 

P, =2^jr — r*— ^"« 
positive or zero. 

7. In forming the derived equation F{x) the only difficulty 
is to form the last term P^. For when this is known the other 
terms can he at once derived/rom it by an easy process. 

Let a, J, c... be the roots of f{x) == with their signs changed, 
and let 

Let A stand for the operation 

da db do 
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Then since — = dbc ... we have obviously 

In the same way we have 
And generally 



P. P. 



and so on up to 



Pii Po ' 



Po 



Let us now operate with A on any expression 

which has the same number of factors in every term. Let r be 
the number of factors, then <f> may be written 

Let P^ P^...P^ be the coefficients of the derived equation 
F(x), and let i^ = 1. Then since the roots of F(x) are the sums 
of the roots oij(x) taken two and two, it is easy to see that 

&c. = &c. 

Thus when P^ is known, the other terms may be calculated 
without difficulty. The term P^ will be called the fundamental 
term of the equation. 

Example. Given in the case of a biquadratic 

P.^PiPJ[>^-PoPz''PiP^ 

to calculate P,. 
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Performing the operation 

Cy e^ Cy t^ Cy m Co 

on Pg we find after division by 2, 

which is the result already given. 

8. It should be noticed that in the equation 

f{x) ^p^ +Piix^ + ... +i?^ = 0, 

if we regard a? as a number, p^p^^ ••• p« are all of equal dimensions. 
It follows from the theory of dimensions, that if any subject of 
operation be the sum of a number of terms of the form 

there must be the same number of factors in every term. For 
example, in every term of the expression for P„ we have 
a + /3 + 7 &c. the same. 

On the other hand, we may regard a; as a quantity of one 
dimension, and in this case Po,p^ '"Pn ^^^^ their dimensions in- 
dicated by their suflBxes. We must therefore have ^H- 27+38+... 
as well as a + ^ + 7+ ... the same in every term. 

These two tests of the correctness of our processes will be 
found convenient. 

9. The whole derived equation being Jcnoion when the funda- 
rn-ental term is knotvn, it is required to find the fundamental term, 

[First Method.] 

If we write — x for a? in any equation, we have a second equation 
whose roots are equal and opposite to those of the first equation. 
If we eliminate x between these two, we shall get a result which 
must be zero when the two equations have a common root. The 
eliminant must therefore contain as a factor the product of the 
sums of the roots of the given equation taken two and two. 

It will afterwards be shown that the last term P^ of the 
derived equation (when p^ is put equal to unity) always contains 
the term 

Pv Pi^ Pi'*' 

with a coefficient which is positive and equal to unity. 
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Hence we have this rule, to find P^, eliminate x between 

oT +jPjaj**"^+jpX'* + ... = 0) 

jPi^""^ +2>8 ^"'' + • • • = ^r 

and divide the result by the coefficient o{pj)j>^..,p^_^. 

It is obvious that the result may be written down as a 
determinant. On trial, however, it will be found more convenient 
to make the elimination by the method of eliminating the highe^it 
and lowest terms than to expand the determinant. 

10. Given the fundamental term of the equation derived from 

f{x)=p^'-'+p,7r'+... + p,.,=o (1). 

to find the fundamental term of the equation derived from 

;?X+;>i^*''+-.+JP«.ia?H-i>n = (2). 

[Second Method.] 

Let Q^ be the product of the sums, two and two, of the roots 
of the equation (2) taken with their signs changed, so that Q^ is 
the same as the fundamental term of the derived equation and 
differs from P^ only in having a suffix more convenient for our 
present purpose. 

Let Q^ be expanded in a series of powers o{ p^ : thus 

where <^o, <^j, <^j, &c., are all functions otp^^ »j, &c., which functions 
have to be found. Let the roots of /(a?) = with their signs 
changed be a, b, c, . ., then 

^»_i=(a + 6)(a + c)... 

Let us introduce a new root, which, when its sign is changed, 
we shall call r. Then 

C„=(a + i)(a + c) ... {r+a)(r + b)... 

jTO 

This value of ^„ must be the same as that given by the series 
when we write 

respectively, for jp^^, p^, p^, &c., p„_^, p^. Equating the coefficients 
of the terms independent of r we have 



^0 = Qn-x „ 



R. A. 
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EquatiDg the terms containing the first power of r we have 

Substituting for ^^ we have, by a known theorem in the Dififerential 
Calculus, 

Equating the terms containing the second power of r we have 

and so on. Thus we have 

If we examine this process, we see that when Q„^^ is known, 
we may at once write down the terms independent of p^ and the 
coefficient of p^. The process to find the coefficient of p^ is 
longer, but it may be much shortened by the consideration that 
when Pq — 1 the result must be an integral function of the coeffi- 
cients. We may therefore omit all terms as soon as they make 
their appearance, which do not contain the factor pVi> ^^^ ^® 
know that such terms must disappear from the result. 

This method is not so convenient as the one which will be 
presently given to find the coefficients of the higher powers of p^. 
But it is useful as showing that Q^ contains the term 

with a positive integral coefficient equal to unity. This will be 
clear from the consideration that the term independent of p^ in 

Q^ is obtained from Q^.^ by multiplying by ^^^ . If therefore Q^^ 

Po 

contains ^'^'';j ^', Q. must contain the term ^^^«';:f "-i . No 

other term can be formed which is equal to this with an opposite 
sign, for the terms which enter by the other processes to be per- 
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formed on Q„_^ all contain |)^ as a factor. Now Q^ = -; therefore 

Po 

0, contains the tenn^^*, ft contains ^'^•A, and so on. 

" Po Po 

It has been shown that every term in p^^Q^ has the same 
number of factors (Art. 8). It follows from this reasoning that 
this number of factors is n— 1. 

11. To find the fundamental term of the derived equation hy 
means of a differential equation. 

[Third Method.] 
The fundamental term required is a factor of the eliminant of 

Let aj' = y, then we have 

py '+;»./*'•+ =o) 

n— 1 t»-3 \ 

i'.r+i'.y^^H- =0 Lodd. 

n-\ w— 3 1 

p^y" +Psy" + =^' 

If we write t/ + dt/ for y the result of the elimination must be 
the same. Hence if we make 






^ n even, 



n odd, 



rfy 2 ^^' dy" 2 P''^''\ 

and if jFbe the eliminant, wo have 

It follows that whether n be even or odd, JFmust satisfy the 
equation 

2—2 
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dE dE 



^^ / dE ^ dE\ 



1 = 0. 



We may make the elimination by multiplying the two equa- 
tions by y, y'..., until we have as many equations as we have 
powers to eliminate. If in the determinant thus formed, we 
multiply out the terms in the diagonal joining the right-hand top 



"-1 ^+1 



corner to the left bottom comer, we get when n is even p * p^ 

n-l n-1 ** 

and when n is odd p^^ * J?o ' • -'^^^ Qn i^^st contain n — 1 factors 

w — 1 

and be of the n — ^ — th degree. Hence when n is even E=- cp^Q^ 

and when n is odd E=cQ^ where c is some constant. 

Now -J— does not occur in the above differential equation. 

Hence treatmg Pj as a constant, we see that Q^ must satisfy the 
differential equation 

V^- dp, + P"-' dp J + [P'-* dp^, + ^- dp J 



+ 3(,„.^^ + ,^-^-.&c. = 0. 



12. We may show that jj^**'* Qn ^^ ^ symmetrical function of the 
coefficients p^, p^...^ p^ and the same coefficients read backwards. 
Let a, 6, c... be the roots of y(a;) =0 with their signs changed, 
then <?« = (« + h) (a + c).... 

If now we read the coefficients in the opposite order, the roots 
of the equation thus formed will, when their signs are changed, be 

-,=-.... If Qn he the fundamental term of the equation derived 
a o 

from this, we have 



^-s-se^)-- 



Since " = ahc... we see that 

prQn=pr'Qn- 
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We may therefore infer that ^o*"*Q^ also satisfies the differential 
equation 

/ dE^ dE\^^( dE ^ dE\^^f dE , dE\ . 

v^d^^p^d^r\p^^^^ 

13. We may use either of these differential equations to find 
Qn when Q^^ is given. 

Let the first differential equation be represented by 

and let 

where A^, A^, ... are functions ofp^, p^, &c. The value of A^ has 
been proved in Art. 10 to be 

Po 

To find the other coeflScients of the powers of p^ substitute 
this value of Q^ in the differential equation ; we have 

+ p^^^A^ + 2p^_^p^A^ 
Equating the several powers ot p^ to zero, we find 

jr«-a 
&C. = &C. 

Thus by one regular and easy process each term may be 
derived from the other. 

In performing this process we may omit every term in the 
subject of operation which does not contain p^j. For p^^ can be 

introduced only by performing -r- , and since p„ is absent from the 

coefficients, this operation yields nothing. 
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lu this way we find 
PoQt =Pi> 

PoQ,'=PiPt-PoP»> 

PoQ, = PiPtP, -PoPs -PiP*' 

P» Qi '=PlP»P,Pi - PtPtP* - PiP4 

-Pi i-P,PJ>» +PiP' - ^PJ>xP^ 

To illustrate this process, consider how Q^ is obtained from Q^. 
The first line is formed by multiplying the line above by p^, this 
is Aq. To find the coeflScient of —p^ we operate with 

on such of the terms in the line above as contain p^ and then 
divide by p^. Performing the same operation on the coeflScient of 

(■"i^ft) we obviously obtain the coeflScient of ^^ . 

In M. Serret's Cours d^Algkbre Sup^rieure^'Note III., there will 
be found a method of forming the last term of the equation to the 
squares of the diflferences, which suggested the method used in 
Art. 13, of substituting in a diflferential equation, if only a diflFer- 
ential equation could be found, [See also Dr Salmon's Higher 
Algebra, Arts. 60, 64«, and 72.] 
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Statement of the tlieory by which the conditions of atahility of a dyna/tnical 
system with n co-ordinates a/re made to depend on 2q conditions. 
Arts, 1 — 4. 

A rule by which these conditions may be derived one from another, to- 
gether with certain other true but not independent conditions. Arts, 

A rule by which, when tJie coefficients of the dynamical equation are letters, 
the 2n conditions of stability may be inferred, one from another, 
wii/iout writing down any other conditions. Arts. 9, 10. 

A metliod by which certain extraneous factors may be discovered and 
omitted. Arts. 11, 12. 

Consideration of the reserved case in which the dynamical equation has 
equal and opposite roots. Arts. 13 — 18. 

A few geometrical illustrations not necessary to the argument. Arts. 
. 19—25. 

Application to a Dynamical Problem. Art, 26. 

1. It has been shown in the first Chapter that the stability of 
a dynamical system with n co-ordinates oscillating about a state of 
steady motion depends on the nature of the roots of a certain 
equation of the 2n^ degree which we may call 

The system is stable if the real roots and the real parts of the 
imaginary roots are all negative. Now Cauchy has given the fol- 
lowing theorem of which we shall make some use. 



Let z^x + y^— 1 be any root, and let us regard os and y as 
co-ordinates of a point referred to rectangular axes. Substitute 
for z and let 

Let any point whose co-ordinates are such that P and Q both 
vanish be called a radical point. Describe any contour, and let 
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a point move round this contour in the positive direction and 

p 
notice how often -^ passes through the value zero and changes its 

sign. Suppose it changes a times from 4 to — and ^ times from 
— to +. Then Cauchy asserts that the number of radical points 
within the contour is ^ (a — ^8). It is however necessary that no 
radical point should lie on the contour. 

2. Let us choose as our contour the infinite semicircle which 
bounds space on the positive side of the axis of y. Let us first 
travel from y = — cx) toy = +a) along the circumference. 



If f{z) =p,z''-\-p/''' + ... -\-p^ 

we have changing to polar co-ordinates 

f{z) =py (cos 71^ + sin w5 V— 1) + . . . 

Hence 

P=py cos nd + py^^ cos (w-1) 0+... ' 
Q = JPo^** sin nO + p/"^ sin (n - 1) ^ + . . . j ' 

In the limit, since r is infinite, 

Tj = cot 710 ; 

P IT 

Y^ vanishes when w^= (2/c + 1) — , i,e. 

^=±-J. ±-l, ±-J (A); 

j^ is infinite when n^ = 2/c ^ , i.e, 

^ = 0, ±?J. ±*J, ±^| (B). 

n z n2 nz ^ ^ 

The values of in series (B) it will be noticed separate those 
in series (A). 

P . 

When is small and very little greater than zero, ^ is posi- 

tive, and therefore changes sign from + to — at every one of the 
values of in series (A). If n be even there will be n changes 
of sign. If n be odd there will be w — 1 changes excluding 

0= ±^, in this case ^ is positive when ^ is a httle less than — , 

TT 

tod negative when ^ is a little greater than ^ . 



and 
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Let us now travel along the axis of y still in the positive di- 
rection, viz. from y= + oo toy = — oo. Since a? = it will be 
more convenient to use Cartesian co-ordinates, we have, since 

and z — ys/—iy 

The condition that there should be no radical point within the 

contour is that this expression should change sign through zero 

from — to + as often as it before changed sign from + to — on 

travelling round the semicircle. If n be even the numerator has 

one more term than the denominator, and when p^ and p^ have 

P . . 

the same sign, j^ begins when y is very great by being negative. 

In order that it should change sign through zero n times, it is 
necessary and sufficient that both the equations 

Pn -Pn^f +Pn^* - . . . = 0, 

p«-iy - Pn^zy^ +Pf^if - . . . = 0, 

should have their roots real, and that the roots of the latter should 
separate the roots of the former. 

If n be odd, the numerator and denominator have the same 

P 

number of terms, and when p^ and p^ have the same sign, ^ begins 

when y is very great by being positive. In order that it should 
change sign through zero from — to + n — 1 times, it is necessary 
and sufficient that the same two equations as before should 
have their roots real, and that the roots of the former should 
separate the roots of the latter. 

In order then to express the necessary and sufficient conditions, 
that f{z) = may have no radical point on the positive side of 

the axis of y, put z = y V— 1 and equate to zero separately the 
real and imaginary parts. Of the two equations thus formed, the 
roots of the one of lower dimensions must separate the roots of the 
other. It is also necessary that the coefficients of the two highest 
powers of z in f (z) should have the same sign, 

3. It has been stated that p^ and p^ the coefficients of the two 
highest powers in f{z) must have the same sign. It is easy to see 
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P 

that, if they had opposite signs, ^ would change sign through zero 

2n times as we travel round the contour. All the radical points 
of the equation would then lie on the positive instead of the 
negative side of the axis of y. 

It has also been assumed that no radical point lies on the 
contour. It has therefore been assumed that/(«) = has no root 

of the form z= 1/ v — 1. It will be more convenient to consider 
this exception a little further on. 

4. It is required to eoopress in an analytical form the conditions 
that the roots of an equation fg (x) = may be all real, and may 
separate the roots of another equation f^ (x) = q/" one degree higher 
dimensions. 

To effect this, let us use Sturm's theorem reversed. Perform 
the process of finding the greatest common measure of /^ {x) and 
j^ {x)y changing the sign of each remainder as it is obtained. Let 
the series of modified remainders thus obtained ^^Q fi^c), f^{x)y 
&c. Then it may be shown that when any one of these functions 
vanishes, the two on each side have opposite signs. It is also clear 
that no two successive functions can vanish unless f^ (x) and f {x) 
have a common factor. This exception will be considered pre- 
sently. 

Hence in passing from a? = — 00 to + 00 no variation of sign 
can be lost except when f^ (x) vanishes. If a variation is lost it is 
regained when x has the next greatest value which makes f (x) 
vanish unless f{x) = has a root between these two successive 
roots of j^ (a?) = 0. Hence this rule: — 

The roots of the equations {^(x) = 0, f^fx) = 0, will be all real 
and the roots of the latter will separate those of the former, if in 
the series 

as many variations of sign are lost in passing from x = — x to 
X = + cx) as there are units in the degree of the equati(m f^ (x) = 0. 

We have supposed the variations of sign to be lost instead of 
gained in passing from a; = — cx) to + x . That this may be the 
case the signs of the highest powers of f^ (x) and f^ (x) must be 
the same. 

These functions are alternately of an even and odd degree, the 
condition that the whole number of variations of sign may be lost 
in passing from a? = — 00 to -f 00 may be more conveniently ex- 
pressed thus : — Hie coefficients of the highest powers of x in the 
series 

f,(x), f,(x). f.(x) ... 

mu^t all have the same sign. 
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6. The process of finding the greatest common measure of 
two algebraic expressions is usually rather long. We may in our 
case shorten it materially by omitting the quotients and perform- 
' ing the division in the following manner. Let 

ft (^) =i>o^** -p,^" + i'X'' - . . . 

/. (^) =i^i^"'' -2>8^"" +P,^"'' - . . . 
then, since p^ is positive, it easily follows by division that 

/3 ix) = ^o;"-' - ^ V-* + ^' V^ - . . . 
where ^^PiP^-PoP^> 

^'='PlP^''PoP6> 

&c. = &c., 

so that by remembering this simple cross-multiplication we may 
write down the value q/" fg (x) witJiout any other process than what 
may be performed by simple inspection. In the same way ^(a?), 
&c. may all be written down. 

6. Ex. 1. Express the conditions that the real roots and 
real parts of the imaginary roots of the cubic 

x^ +px^ -i-qx + r^O 
may be all negative. 

fi{x) = x^-qx, 

/,(a;)=_pjj'-r, 

Mx)=(pq-r)r. 

The necessary conditions are that 

p, pq-^-r and r 
must all be positive. 

Ex. 2. Express the corresponding conditions for the bi- 
quadratic 

x* +px^ + qx* -{-rx + s = 0, 

f{oc)=x^ -qx'^ + s, 
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The conditions are that 

p, pq — r, {pq — r) r-'p^s and s 
must be all positive. 

These are evidently equivalent to the five conditions that 

P, q, r, 8, {pq-r)r-p% 
should be all positive. 

In both these examples all the numerical work has been 
exhibited. 

7. Since the coefficients of the highest powers of x in/Jx) and 
f^{x) Sire Pq and jp^ we see that the condition that^p^ and p^ should 
have the same sign is included in the general statement that all 
the coefficients of the highest powers should have the same sign. 
If the function /(a?) be of n dimensions we thus obtain n necessary 
and sufficient conditions. 

On examining these conditions in the cases of the cubic and 
biquadratic it will be seen that they cannot be satisfied if any one 
of the coefficients of the given equation should be negative. 

8. Although the theorem in its present form gives n conditions 
as the proper number for an equation of the n*^ degree, yet it is 
important to notice that it gives other conditions also which are 
true and may be useful. It has been shown in the second Chapter 
that all the coefficients of the equation f(x) = must be positive, 
hence the roots of^(a;) = must all be positive. It may be shown 
also, that the roots of each of the functions f,{x), f^(x), &c. are 
separated by the roots of the function next below it in order. 
Hence the roots of all these functions must be positive, and there- 
fore in every one of the functions the coefficients of all the powers 
must be alternately positive and negative and not one can 
vanish. If however f^(x) and f^{x) have one or more common 
factors some of the {unctions fjx),fj(x), &c. will wholly vanish. 

9. When the degree of the equation is very considerable 
there is some labour in the application of the rule given in Art. 5. 
The objection is that we only want the terms in the first column, 
and to obtain these we have to write down all the other columns. 
We shall now investigate a method of obtaining each term in the 

first column from the one above it without the necessity of writing 
down any expression except the one required. 

We notice that each function is obtained from the one above it 
by the same process. Now 

/a(^) = (PiP« -AP3) ^"^ - (;>iP4 -PoP^^' + . . . 
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The first and second lines will be changed into the second and 
third lines by writing for 

P^ P^ &c. 



Po. 


P.. 


the values 




Pv 


PiP.-PoP.. 



If then in any term of any function we make these changes, we 
obtain the corresponding term of the function next in order. 

10. Example. Express the conditions of stability for the 
quintic 

f(x) ^p^a? ■\-p^x'' -^-p^x^ •\-p^x^ + P.^ +l>a- 

We have 

/i(«)=2>o^+--- 
/8(^) = (PJ>2-PaP8)«'» 

/4(«^) = {(PlJP2-AP>8-Pl(i>lP4-PoP5)K' 

/sC^) = i[{pj>2 - Poi>8);>8 - Pi (PiPa - PoPb)} (PiPi - P0P5) 

-iPiP^-PoP^Y-P^^* 

fs(^)='[{(PlP2''PoP^Ps-PiiPlP^-PoP5)](PlP4-PoP,) 

- {P1P2 - PoP^^ .i>J {PJ>^ -PoP3)P,' 

11. On examining the conditions as given in the cases of a 
biquadratic and quintic it will be apparent that several contain 
the previous conditions as factors. Thus the analytical expressions 
are rendered much longer than is necessary. It is now proposed 
to investigate a method of discovering and omitting these extraneous 

factors 05 they occur, and thus obtaining the required conditions in 
their simplest forms. 

Let the coeflScients of the several powers of x in the functions 
be when taken positively 



fli^)=p<^ 


Pi' 


Pt' 


Ps'- 


/»(^)=P1. 


P» 


P» 


p,... 


A{x) = A, 


A', 


A", 


A"' 

^x • • • 


fM^B, 


B', 


B", 


■Dill 
JD ... 



&c. = &c. 
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Let US first find which of these terms contain pi as a factor. 
Putting ^1 = and using the rule in Art. 5, the series become 



Po. 


r». 


^4. 


p» 


0, 


Pz> 


^M 


p,. - ■ 


■PoPaJ 


-PoPs' 


-PoPt 


-P(,p» ••■ • 


PoPt> 


-PoPiPt' 


~PoPaP7' 


-PoPsPo 


0, 


0. 


0. 


0, ... . 


0, 


0, 


0, 

&c. 


0. ... . 



Hence the C's and D's all vanish and therefore contain jp^ as 
a factor. By the rule in Art. 5, the E'a contain p^, the Fs con- 
tain Pj", the G'ap^, the ITap^, and so on. 

But since each line is formed from the preceding by a uniform 
rule, it follows that the D*s and jE's contain ^ as a factor, the 
F'a contain A\ the Q*s contain A^, the jET's contain A\ and so on. 

The factor A in the D's and E*a has its origin in the factor p^ 
which occurs in the (7's and D's and would not appear if that 
factor had been omitted when the (7's and i)*s were formed. The 
factor A in the D's and jBTs in the same way gives rise to the 
factor B in the jE's and F's. So that if we take care each time 
we perform the process described in Art. 9 to omit the common 
factor p^ whenever it occurs, all these subsequent factors will never 
make their appearance. 

We shall now show that if these factors are omitted, the 

w — 1 

dimensions of the w** function ^(ar) will be n —^ — . First consider 

the actual dimensions of each function before the factors are 
omitted. If we examine the rule by which each function is 
derived from the preceding, it will become evident that, the 
dimensions of each letter being indicated by its suflSx, the 
dimensions of any function are equal to the sum of the two pre- 
ceding + 2. 

In the following table the first column indicates the function. 

In the second column will be found the dimensions of the leading 

coeflScient of that function when calculated by the rule in Art. 6. 

In the third column will be found the dimensions as given by the 

n — 1 
formula n — ^ . In the remaining columns are the dimensions 

of the extraneous factors ^j, Ay J5, &c. introduced into each term. 
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fM 














fM 


1 


1 








fM 


3 


3 








M<^) 


6 


6 








fM 


11 


10 


1 






.W) 


19 


15 


1 3 






/, (^) 


32 


21 


2 3 6 






/«(^) 


53 


28 


3 6 6 


10 




/,W 


87 


36 


5 9 12 


10 


15 


&c. 


&c. 


&c. 


&c. &c. &c. 


&c. 


&c. 



Each term in the second column is the sum of the two terms 
just above it + 2. The w*^ term in the third column is equal to 
the term just above it + (ti - 1). In all the other columns each 
term is the sum of the two terms just above it. The last term in 
the n*^ row is equal to the (n—S)^ term in the third column. 
We wish to show that any term in the second column is equal to 
the sum of all the terms in the same row to the right of that term. 
It is not diflScult to show from the data just given that if this be 
true for any two adjacent rows, it is true for all the others, and 
hence we may assume it to be always true. 

It is clear that these extraneous factors may be omitted since 
by the conditions already expressed they are all positive. When 
omitted as they occur, the dimensions of the n^^ function has just 

been shown to be n — »— . It is easy to see that the conditions 

thus reduced must contain the terms 

Pi* PiPv PxVJPv TxViVzVi," ^^^ 
Now if we take any one of these as 

PJ>J>zPiP^Pf^ (1)> 

and operate by the rule in Art. 9, we have 

(PiP^ - PiP^P^ (PiPi - PoP5)P, (PiPe - PoP,)Pv 
which contains the term 



PiPiPfJPiPiP^PxPiPi (2)- 

Thus we have p^ introduced as often as there is a factor in (1) 
with an odd suflSx. But it should be introduced only once. These 
extra ^^'s are the extraneous factors to be omitted. Each of these, 
if left, would appear as the factor p^p^ —PoPs ^^ *^^ ^^^^ condition, 
and be still more complicated in the next after that. 
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In order then to obtain the several conditions in their simplest 
form it is only necessary after performing the operation described in 
Art. 5 or Art. 9 to divide by p^^ where k is one less than the number 
of factors with odd suffixes in the condition operated on. 

12. Example. Express the conditions of stability for the 
sextic 

f(x)=p,x'+p^a^'{-p^x'+p^x'' + p,x''-^p,x + p,^0. 

We have 
f(x)=p^'' + ... 
f{x)=p^x'+... 

f {^) = {P1P2P3 - PoPs - Pi Pi +PoPiP,) aJ* + . . . 

/s (aj) = {PlPJ>^Pi - PoP.'Pi -PlPi + ^PoPiPiP5 •-'P1P2P5 +P0P2PSP6 

-PoP^ +l>i>2P6 -PoPiPf^Pe] 35* + ... 

A (^) = [PlP2P3PiP5-PoP,^PiPf^-Pt*PiP6+ ^PoPiPiP^-PlP^P,^ 

+ PoPJ>J>,* - i>oV + ^PiPJ>iP6 - SPoPiPsP6P6 -PiP^PsPe 

^PoPsPe+PiPsPiPe -PiPe] ^ + &C. 

fi (^) ~ coefficient of a; in ^ a; x by p^. 

13. In the preceding theory two reservations have been made. 

1. In applying Cauchy's theorem it has been assumed that 
there were no radical points on the axis of y. 

2. It has been assumed that P and Q have no common 
factor, so that none of the functions/,, /,,, &c. vanish absolutely. 

If any radical point lie on the axis of y, it is clear ih8itf{z) = 
must have a factor of the form (z^ + a')*". Let f{z) = {z* + a'*)*' <l>{z). 

In this case when we put ;3 = y V — 1, we have 

Thus P and Q have a common factor, and we are warned of the 
possible existence of radical points on the contour by the total 
vanishing of some one of the functions /,, j^, /, &c. 

The two reserved cases may therefore be included in one. If 
f{z) = be the equation furnished by dynamical considerations, 
we form the functions/,, /„ /„ &c. If all these be finite, the 
question of the stability of the system has been answered. If any 
one vanish absolutely, / (a?) and / (x) have a common measure. 
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and we must add some further considerations. It will be con- 
venient to examine separately the dynamical effect of the roots 
which do not and which do enter through the greatest common 
measure. Let us begin with the former. 

14. Following the same notation as before, we have 

P^±My)=Pn-Pn^+ 1 

If then / (z) have two roots, viz. ± (A + A; V — 1), which are equal 
and opposite, then f^ (y) and /, {y) must have two common roots, 

viz. + — . The common measure therefore of f^ (y) and 

f% (y) contains all the roots of /(y V — 1) which are equal and 
opposite. Conversely the greatest common measure of P and Q is 
necessarily an even function of y, and if it be equated to zero, its 
roots are necessarily equal and opposite. These roots must also 

satisfy/(yV^) = 0. 

Let this greatest common measure be yjr (y^) = 0, and let j^ be 
the highest power which enters into it. Also let 

then ^ (-a^) is a function which, as has just been shown, has not got 
two roots equal and opposite, and to this function we may apply 
Cauchy's theorem without fear of failure. Putting « = y V — 1, let 

P' 

Then we wish to express the condition that ^y should change 

sign from — to + through zero w — 2/^ times if n be even and 
u — 2r — 1 times if w be odd. But 



/(yV-i) = P+QV:ii, 

and /(y>/"=n[) = t(2^(^ + e'V^); 

/. P = ^(y»)P'| 

Thus the number of changes of sign in -^ is exactly the same 

p 

as that of -q . The factor -^ (y*) will run through all the functions 

/s (y)» fs (y)> ^^- obtained from f^ (y) by a process which is equiva- 
lent to that of finding the greatest common measure of /^ (y) and 

B. A. 3 
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/ (y). The changes of sign of this factor will therefore not affect 
tne number of variations of sign in the series /i> /j, /g* &c. 

The last factor which is not zero is/^j.^^ (y) if n be the dimen- 
sions of/ (y). 

'Rence if we omit the considerations of the vanishing factors and 
apply the same rule as before to the 11 + 1 — 21 remmning factors, we 
can escpress the condition that the proper number of changes of sign 
from — to + have been lost through zero in the function <f> (z),i.e. that 
the roots not given by the vanishing of fn+j^ are aU of the character 
to ensure stahility, 

15. Let us next consider the effect on stability of the roots 
indicated by the absolute vanishing of one of the subsidiary func- 
tions. This function must be of the form 

where n is even. The corresponding factor of f{z) is 

The roots of this equation are two and two equal with opposite 
signs, it is therefore necessary for stability that no root should 
have any real part. To express this condition, draw a straight 
line parallel to the axis of y at an indefinitely short distance from 
it, viz. a? = A. Let us apply, in the same manner as before, Cauchy's 
theorem to the contour formed by this straight line and the infinite 
semicircle on its positive side. Putting 2^ = A + «', we have 

F{z) = y. + 2q^Jiz' + q^ji'' + ^q^hz'' + ... 
The two functions are therefore, omitting the positive factor h, 

Now/ (y) and/j, (y) are what P and Q become when arranged 
in descending powers of y and the coefficients of their highest 
powers Toade to have the same sign. Hence 

The rule described in Art. 4 will now become the same as 
that usually called Sturm^s theorem. We are to seek the greatest 
common measure of/ (y) and its differential coefficient, and make 
the coefficients of the highest powers of y in these two and in the 
series of modified remainders all positive. 
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That we should have been led to Sturm's theorem in this case 
is just what we might have expected. For to express the condi- 
tions that the roots of 

are all of the form ± r V — 1 is the same thing as to express the 
conditions that the roots of 

are all real. 

16. There is however another mode of proceeding. Suppose 
we have calculated the functions f^yf^y &c. for the general equa- 
tion 

and find when the values of p^,, jj^, &c. are substituted that some 
one function say f^ of the series absolutely vanishes, and therefore 
also all the functions which follow it. Then operate on each of 
these vanishing functions with 

repeating the operation until we obtain a result which is not zero. 
If we now replace these vanishing functions by these results we 
may apply the rule of Art. 4, just as if these were the functions 
supplied by the process of the greatest common measure. As this 
process is not so convenient as that already given it is unnecessary 
to consider it in detail*. 

17. As a numerical example, let us examine whether the roots 
of 

/ (a?) = a^ + 2a;' + 4aj* + 4^' + 6aj* + Ga' + 7a;* + 4a; + 2 = 

satisfy the conditions of stability. In order to show the working 
of the method it will be necessary to exhibit all the numerical 
calculations. We have by Art. 5, 

* The fanotion /^ (a;) vanishes beoanse the equation f{z)=0 has two roots equal 
and opposite. If we put z=i^-^hj where ^ is as small as we please, this peculiarity 

will disappear. Thus if the values of z are of the form ±(oi/3 J-l) tiie corre- 
sponding values of ;;' are -h±a±p ^J -1. These values of z wiU indicate stability 
if a be zero and instability if a have any value positive or negative. If ^ be as smaU 
as we please and positive, the v^ues of z' wiU indicate stability or instability under 
the same circumstances. We may therefore apply the rule of Art. 4 to the func- 
tion f{s^+h) instead of /(2), provided we retain only the lowest powers of h which 
occur. Hence all the functions fi{x), /•(«)... /r-i(«) which do not vanish are 
unaltered. To find what function wiU replace fr {x) we must increase by ^ aU the 
roots of /(z)=0 when their signs have been changed. This may be effected by 
performing on / (z) tiie operation represented by A in Art. 7. The rule in the text 
therefore follows from the one given in that article. 

3—2 
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/,(ir) = 2aj'- iaj' + eaj*- 4», 

/ {x) = 4a:^ - 4a^ + 8a?, 
/. (a?) = 8«* - Sa:* + 16. 

Here we find/g (a:) to be absolutely zero, accordingly by Art. 15 
we replace it by the diflferential coefficient of /^ (a?), this being 
Sturm's rule. We have therefore 

/,(aj)=8(4a^-24 

f,{x)=^^8\20.x, 
};(a?) = -8M20. 

We see that the two last of the coefficients of the highest 
powers are negative. The roots therefore do not satisfy the con-^ 
ditions of stability. 

As another example, take the equation 

/a: = a:« + a^ + 6a?* + 5a? + 1 la:* + 6a? -I- 6. 

Here 

/^(aj) = a;'-6a;*+ll«"-6, 

/, (a?) = a?" - 5a?» + 6a?, 

/,(a:) = a?*-5a?^+6, 

/,(a?) = 0. 

Keplacing/4 (a?) by the differential coefficient of ^(a?), we have 

/4(a?) = 4a?^-10a?, 
f, (a?) = lOa?^ - 24a?, 
/, (a?) = 4a?, 
/(a?) =4.24. 

Here all the coefficients of the highest powers are positive, 
hence the roots satisfy the conditions of stability. 

It is clear that when the coefficients are numerical the rule 
given in Art. 5 is the most convenient, but when the coefficients 
are letters, the rule in Art. 9 will be found preferable. 

The process would be simplified by omitting the alternate posi* 
tive and negative signs of the terms in each line. 
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18. It' may be interesting to express the two subsidiary funo- 
tions fj (x) and f, (x) in terms of the roots of the given equation. 

Let ttj, a,, a,...a^ be the roots of the given equation/(a?) = 0, so 
that 

/(a?) = (a?-aO [x-a^ipc-a^... 

= af" + p^af'^ + p^x"^ + ... 
Then it is evident that 
±/j {x V^) = aj* +p^'^ +Pi^'* + . . • 

= J (x + a^) (x + a^) ... + H^"-«i) (^-^s) •- 
It may be shown* that 

19. The foUowing propositions are not necessary to the main 
argument, but as they illustrate geometrically the propositions in 
this chapter it has been considered proper to state them very 
briefly. The demonstrations wiU therefore be much curtaUed. 

The equation being 

we put as in Art. 1, z = x + y V— 1, and thus obtain two curves, 
whose equations expressed in polar co-ordinates are 

P = j[?^r* cos 710 +py-^ cos (n - 1) ^ + . . . = 0) 
g=;)^r*sinw^+i?/*"'sin(7i-.l)^ + ... = 0J • 

These intersect in the radical points of the equation/ (2^) = 0. 

20. If we trace these curves we find that the curve P = 
has n asymptotes whose directions are given by cos n0 = 0, ie, 

>j__l7r Stt 5 it 

w2' w2' n2' 



* Let ns assnme 

where A^, A^^ &c, are constants whose values have to be fonnd. Patting 0;=%, 
we have jPi%*~^ +JP2^i*~' +&c. = Ai (a^ - a^) . . . {ai - a J. 

But since a^+jPia{»~^ + ... = (a;-ai) (x-a^)...{x-aj, 

we have by patting x^a^ and x= -a^^ 

Oj* -l>iai*~^ + .. . = 2ai («! + ttg) . . . (Oj + a^). 

Subtracting the- second of these results from the first we find A^to have the value 
given in the text. 
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These asymptotes all pass through the same point on the axis 

of X, viz. a; = — -^ . It is also clear that only one branch of the 

curve can go to each end of an asymptote. Similar remarks apply 
to the curve ^ = 0, the directions of its asymptotes being given by 
sinw^ = 0, i.e, 

/}_f\ 2 TT 4 TT 6 TT 

^"" ' n2* n2' w2' 

21. From these simple propositions we might, if it were 
worth while, deduce that every equation must have a root. The 
asymptotes of the two curves P = 0, Q = are alternate, and no 
two branches of the same curve can approach the same end of an 
asymptote. By sketching a figure, it may be easily shown that 
some branch of the P curve must cut some branch of the Q curve. 

22. Let us next consider the intersections of the curves P = 0, 
Q = 0. 

If we transform the origin to h, k, we put x='h + ^,y=^k + rj. 

This is the same as expanding f{h + A;V— 1+f + ^V— 1), and 
collecting into two parcels the real and imaginary terms. Let the 
expansion be 

where A^A^... are of the form 

c (cos a + sin a V^^). 



If we put f + ^ V— 1 = r (cos ^ + sin ^ V— 1), we have 

P = Co cos a^> + c^r cos (^ + aj + c/' cos (2^ + ag) + ...| 
^ = Cq sin a^ + c^r sin {0 + a^) + c^r^ sin (2^ + a J + ... j 

If the point {h, A?) be a point of intersection c^ = 0. If the inter- 
section be a double point on either curve, the terms of the first 
degree must be zero, therefore c^ = 0, and the origin is therefore 
a double point on the other curve also. 

It is not difficult to show that if the intersection be a multiple 
point of any degree of multiplicity on one curve, it is a point of 
the same degree of multiplicity on the other curve. The tangents 
to these branches all make equal angles with each other, the 
tangents to the P and Q curves being alternate as we travel round 
the point of intersection. If the intersection be not a multiple 
point on either curve, the branches cut at right angles. 

Let us travel round a point of intersection along the circum- 
ference of a small circle whose centre is the point of intersection 
in the direction in which is measured. Then it may be shown 
that as we pass from a P curve to a Q curve, P and Q have 
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opposite signs, and aa we pass from a Q curve to a P curve the same 
sign. This is in fact merely Cauchy's rule for the changes of sign 



23. Let US express the condition that there is no radical point 
on the positive side of the axis of y. This is the geometrical 
illustration of Art. 2, 

Draw a circle of infinite radius, and let it cut the asymptotes 
of the F curve in P,, P^P^,..,P^ and the asymptotes of the 
Q curve in Q^, Q,,.--^^' These points alternate with each other. 
Taking only those points which lie on the positive side of the 
axis of y, the P and Q curves may be said to begin at these points 
and are to intersect each other only on the negative side of the 
axis of y. The branches of the two curves must therefore remain 
alternate with each other throughout the space on the positive 
side of the axis of y. Their points of intersection with the axis 
of 1/ must be also alternate, and hence if we put a; = 0, in the 
equations P=0, Q = 0, and regard them as equations to find y 
the roots of each must separate the roots of the other. 

Conversely, we may show that if the intersections of the two 
branches are alternate on the axis of v, they cannot have intersected 
on that side of the axis of y on which the common intersection of 
all the asymptotes is not. This is the result arrived at in Art. 2. 

24. The following diagrams exhibit the forms of the curve 
p=0, Q = for a biquadratic. The dotted lines represent the 
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A biquadntie with eqoftl imsginarj root*. 
26. [If «=a: + y V — i, we have 

Bifferentiatiiig with respect, firstly to x and secondly to y, we 






^+Svz 



^1 



'- and ^ = -^; 30 that hoth 
ay 



da dy dx 

the functions P and Q satisfy the equation 

The equation y^z) = gave ua two curves which we have called 
P = Oand Q=0. In the same way the derived equation/'{s) = 
will give us two other curves, which we may represent by P' = 
and Q' = 0. These we may call the derived i* and Q curves. 

If as in Art 22 we transform the origin to the point {k, k) we 
have 

P' = c,cosa,+2(yco8(e+aJ+3c/co8(2e + a,l+... 
Q' = c, sin 5[, + 2c^ sin {^ + oJ+ 3c,i* sin (25 + o^ + . . . 

If the origin he at a point of intersection of the curves P = 0, 
P = which is not a double point on the first of these curves, we 
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have Co cos Oo = and cos a^ = 0. Hence the tangent to the curve 
P = at the point of intersection is parallel to the axis of a?. 
Conversely if we move the origin to any point on the curve P = 
at which the tangent is parallel to the axis of x, we find that the 
curve P' = passes through the origin. Hence the derived P 
curve passes through all those points on the curve P = at which 
the tangent is parallel to the axis of x, and all those points on the 
curve (3—0 at which the tangent is parallel to the axis of y. In 
the same way the derived Q curve passes through those points on 
the curve Q = at which the tangent is parallel to the axis of x, 
and those points on the curve P=0 at which the tangent is 
parallel to the axis of y. 

If Co cos a^ = 0, and Cj = 0, the origin is a double point on the 
curve P= and the origin also lies on both the curves P" = 0, 
Q =. 0. So that both the derived curves cut the curves P = and 
Q = in their double points. In other words, these double points 
are radical points for the derived equation /'(^a?) =0. In the same 
way, any multiple point on either curve is a multiple point of one 
degree less multiplicity on both the derived curves. 

If a finite straight line AB be drawn parallel to the axis of x 

joining two points Ay B on the same or on different branches of 

the curve P = 0, this finite straight line must cut one or more 

branches of the derived P curve. For it is clear that if P vanishes 

dP 
at A and P, P' which is equal to -7— cannot keep one sign between 

A and P, and must therefore vanish somewhere between A and B. 
If A and B be adjacent points, i,e. if there be no other points 
between A and B belonging to the curve P = 0, then the straight 
line AB must cut an odd number of branches of the derived P 
curve. In the same way if a straight line CD be drawn parallel 
to the axis of y joining two adjacent points on a derived Q curve, 
this straight line must also cut an odd number of branches of a 
derived P curve between G and D. 

By considering a tangent as the limit of secant, it again follows 
that if a tangent be drawn to the curve P = parallel to the axis 
of a?, the derived P curve must pass through the point of contact. 

Let 5 be a radical point on the derived curve /*(«) = 0, and let 
it not be a double point on either of the curves P= 0, Q = 0. Let 
a straight line be drawn from R in any direction cutting the 
branches of either of the curves P= 0, Q = in the points A^^ A^ 
&c. Then we may show that 

so that the polar line of R with regard to either of the curves 
P = 0, Q = is at infinity. 
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The positions of the radical points of the derived equation 

Z' («) = relatively to any branch or branches of the curves P = 0, 

Q = may be found by the use of Cauchy's theorem. If the 

point spoken of in Art. 1 travel along a branch of the curve P=0, 

it is easy to see that 77 = ;t^ • ^ li travel along a branch of the 

curve Q = we have 7y = — t~ • If *ten any contour be partly 

bounded by branches of these curves, the simplest inspection of 
the points at which the tangents are parallel to the axis will 

determine the changes of sign of 7y as it passes through zero. 

If another part of the contour be an arc of a circle of infinite radius 
whose centre is the origin, the changes of sign through zero will 
be from + to — and their number will be indicated by the number 
of asymptotes of the derived P curve which cut the arc] 

26. The use of Watt's Governor in the steam engine is too 
well known to need description. It has however, as commonly 
used, a great defect. It is sometimes of importance that the 
engine should continue to work at the same rate notwithstanding 
great changes in the resistances. Suppose the load suddenly 
diminished, the engine works quicker, the balls diverging cut off 
the steam, and the engine, after a time, again works uniformly, 
but at a different rate from before. The balls as they open out 
or close in are usually made to describe circles. Let them now be 
constrained to describe some other curve which we may afterwards 
choose so as to correct the above defect. If this curve be a parabola 
and the balls be treated as particles, it is clear from very elemen- 
tary considerations that these will be in relative equilibrium only 
whin the engine works at a given mte. This pknciple is dui 
to Huyghens, see Astronomical Notices, December, 1875. It is 
now proposed to determine the condition of stable oscillation about 
a state of steady motion. 

Two equal rods AB, AB are attached at A by hinges to a 
small ring which can slide smoothly along a vertical axis. The 
ring is attached by a rod to the valve and can thus govern the 
amount of steam admitted. Two equal balls are attached at B 
and By and the centre of gravity O of the rod AB and the ball B 
is constrained to describe some curve. To represent the inertia of 
the engine we shall suppose a horizontal fly-wheel attached t6 the 
vertical axis whose moment of inertia about the axis is /. Let the 
excess of the action of the steam over the resistance of the load 
be represented by some couple whose moment about the vertical 
axis is /(^), where Q is the inclination of the rod AO \^ the 
vertical, and / is a function which depends on the construction of 
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the engine. Since the steam is cut off when the balls open out, 
it is clear that -^^^ is negative. 

There may be also some resistances which vary with the 
velocity. Let these be represented by a couple B -^ tending to 
retard the motion round the vertical axis and a couple round A in 
the plane BAB^ equal to mO -^ ^ where ^ is the angle the 

vertical plane BAS makes with a fixed vertical plane. 

Let m be the mass of either sphere and rod; k the radius of 
OTration about an axis throus^h perpendicular to the rod, and k' 
E .bout .be r<A^ I. la ffi «,. equation of .«guta 
momentum gives 

Let the steady motion be given by ^ = a, -~ = w, and let the 
oscillations be represented by ^ ^ a + a?, -^ = n + y. We have then 

The equation then reduces to 

{/ + 2m (;?T? sin* a + P cos* o)) ^ + 2m (? + &• - fc'^ sin 2aw J 

=f{a)x-Btf. 
This equation may be briefly written 

Let z be the altitude of above some fixed horizontal plane. 
Then if T be the semi vis viva 

2 r= J^'* + 2m {(ifc* + Z^ sin*^ + U* cos»5} c^" + 2mA;»^ 

+ 2mj(jy+Pcos»^U'*. 

If U be the force function, omitting the couple of resistance, 
we have U='-2mgz. The virtual moment of the couples of 
resistance being 2mCff&6, the Lagrangian equation of motion 

becomes 

d dT dT dU ^ ^^ 

^^Wdd^Te"^'^^^' 
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Substituting for Tand Uwe have when the motion is steady 

!L'(p + A;«_jk'«)sin^cos5 = ^, 



• • 






Since Z sin ^ is the distance of from the vertical axis, we see 
that the path of G must be a parabola. The semi latus rectum is 

5-^ 



f2\ 9 



which, we notice, is independent of the radius of the balls. The 
length of this latus rectum must of course be adjusted to suit the 
particular rate at which the engine is intended to work. 

When the system is oscillating about the state of steady motion 
we have, putting 

a« = P + ^(Z" + A;*-ifc'ysin*a, 

and rejecting the squares of x and y, 

(JK» + a«cos»a)^ + (7^-(P + A;'-Jfc'«)sm2wy = 0. 

The term x, it will be noticed, has disappeared from the 
equation. This equation may be briefly written in the form 

Eliminating y from the equations of motion we have 

The coeflScients are all positive, the necessary and suflEicient 
condition of stability is therefore 

{AC+ BH){BG+EL)>AFHL. 

In some clocks to which Watt's Governor is applied, there is a 
special arrangement which causes (7 to be much greater than B. 
See the Astronomical Notes, xi., 1851, and the Memoirs of the 
Astronomical Society^ Vol xx. Neglecting therefore B, we have 

CE>FH. 
.-. 2(7mn(P + ft'-ifc'^sin2a>^(A?" + a"cos"a). 
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Formation of the equationa of steady motion and of small oscillation 
where Lagrang^s method may he used. Arts. 1 — 5. 

The equations being aU linear the conditions of stability a/re expressed 
by the character, of the roots of a determina/ntal equation of an even order. 
Art. 6. 

Mode of exfpamding the determinant. Art. 7. 

A method of finding the proper co-ordinates to make tite coefficients of 
the Lagrangia/n function constant. Arts. 8 — 10. 

ffow the Hwtnwnic oscillations about steady motion differ from those 
about a position of equilibrium. The forces which cause the difference are 
of the nature of centrifugal forces produced by an im>aginary rotation 
about a fixed straight line. Arts. 11 — 19. 

Reduction of the fundamtental determinant to one of fewer rows by the 
etvmvnaMon of all co-ordinates which do not appea/r except as differential 
coefficients in the Lagrangian function; vnth am, example. Arts. 20 — 23. 

Formation of the equations of Motion a/nd of the determinant when 
the geometrical equatioTis contain differential coefficients, so that Lagra/rvg^s 
msthod cannot be used; with an exa/mpU. Arts. 24 — 27. 

1. Let the system be referred to any co-ordinates f , 17, ^ &«• 
The general expression for the kinetic energy is 

where P, Q, &c. are known functions of f , 17, f, &c. and accents 
have their usual meaning. Let us suppose the system to have 
some motion represented by . 

l=/(0> V^F{t\ &c. 

and when disturbed, we wish to find the oscillations about this 
motion. To eflfect this, we put 

f=/(0-f^, V = F{t) + <l>, &c. 
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where 0, ^, &c, aire all small quantities. Substituting and expand- 
ing T in powers of 0, ^, &c., we find 

+ ^{AJ' + 2AJ<I>+...) 

+ CJff + CJ<l>'+C^e'<l> + ... 

+ &c. 

In the same way we may make an expansion for the Potential 
Energy of the forces, viz. 

+ ^{E,^0' + 2EJ<l> + ...l 

when these two functions are given the whole dynamical system 
and the forces are known; and we may form the equations of 
motion by Lagrange's method. 

2. We shall here however limit the question by supposing 
that the motion about which the system is oscillating is what has 
been called in Chap. L steady. The analytical peculiarity of such 
a motion is that when referred to proper co-ordmates, every coeffi- 
cient in each of these two series is constant, i.e. independent of t 
As already explained the physical peculiarities are that the vis 
viva is constant throughout the steady motion and the same oscil- 
lations follow from the same disturbance at whatever instant it 
may be applied to the motion. A method of discovering the 
proper co-ordinates, if unknown, wiU be given a little further on. 

3. In order to form the equations of motion we must now 
substitute in Lagrange's equations 

dtdff d0^ d0' ' 

&c. = 0, 

rejecting all the squares of small quantities. The steady motion 
being given by 0, ^, &a all zero, each of these must be satisfied 
when we omit the terms containing ^, ^, &c. We thus obtain the 
equations of steady motion, viz. 

^a = K* 
&c. = &c. 
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These equations may be simply formed in any ease by the 
following rule. Putting L = T—V, so that L is the difference 
between the kinetic and potential energies, expand L in powers of 
the co-ordinates 6, ^, &c. regarding ff , ^', &c. as zero. The required 
relations are obtained by equating the coefficients of the first 
powers to zero. This rule may be also expressed thus. Let L he 
the general expression for the excess of the kinetic energy over the 
potential energy of a dynamical system in terms of its n co-ordinates 
^, 97, &c. Let this system he mxyving in steady motion with constant 

values 0/ jf > ;n > &c.* Then substituting these constant values in 

the general expression for L, the relations hetween the constants of 
steady motion are given hy 

dL ^ dL ;. « 

In this way we obtain in general as many equations as there 
are co-ordinates. Usually the coefficients in the expression for T 
are constant because some of the co-ordinates are constant in the 
state of steady motion, and the other co-ordinates appear in the 
expressions for T and V only as differential coefficients. In such 
cases we have clearly fewer equations than co-ordinates to con- 
nect the constants of steady motion. We have then a system of 
possible steady motions which we may conveniently term parallel 
steady motions. 

4. To obtain the equations to the oscillatory motion, we retain 
the first powers of 0, ^, &c. We thus obtain a series of equations 
of which the following is a specimen : 

+ ...=0. 

To solve these we write 

e^M^e"^, <l>=:M,e'^, ^^jjf.e'^, &c. 

Substituting we obtain on eliminating the ratios M^: M^i M^ &c. 
a determinantal equation, viz. 

* Since we may in Art. 1 change the co-ordinates from ^, 17, I*, &o. to ^], 17^, j'^, 
Ac., where ^=/(f j) 17=2^(171) &o., it is clear that the steady motion can be always 
expressed by constant values of the differential coefficients of the co-ordinates. 



* 
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By^m*-A^+E^ 


B^*-A^ + Ea 


B„m*-Aa+Ea 


&c. 




B^* -A^ + E„ 


B^m'-A^+E^ 


&c. 




B„m*-A^ + E^ 


B„m' -A„ + E^ 


&c. 


&c. 


&c. 


&c. 


&c. 



= 0. 



This equation will be referred to as the Determinantal equa- 
tion. 

5. If we refer to the equation formed by this determinant 
and read it in horizontal lines, we have of course the several 
equations of motion, each term being the coefficient of 0, ^, '^, &c. 
in order. In this form the equations may be reproduced by the 
following easy rule. 

Taking the expression for T- F as given in Art. 1, let us 
consider only the terms of the second order, those of the first 
having been already used to determine the steady motion as ex- 
plained in Art. 3. Separate from the rest, the even powers of 
^, ff> ^, i>\ &c. and write for ^^ ff<f>', &q. - 2>*e", - 1^0^, &c., so 

that D will stand either for -j- ox, for the m in the determinant 

at 

of Art. 4, when we write = M^e*^, <f> = M^e"^, &c. Let the sum 

of these terms be called P, so that 

P = i (A,, - E,, - BJT) ^ + ( J„ - E,, - B^D^ 0<l> + &c. 

Let the remaining portion of the terms of T— F, viz. those con- 
taining both 0, 0, &c. and ff, ^', &c., be called Q, so that 

Q=CJff^CJ4;^C^^<t>0'+... 
Then the several equations may be formed from the rule 

dP,dQ_ndQ_Q 

dff^de ^ dff~^ 

dP_^dQ_j)dQ^Q 
dip d<f> dcj) 

&c. =0 

In applying this rule no accented letters will occur except in 
the second term of each equatioa If we wish D to stand for m 
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in the detenninant, we must regard ^, <f)\ &c. as abbreviations for 
BO, D(f>y &c. If we wish to use the equations themselves, we 
replace 1)0, D<f), &c. by ff, ij>, &c. 

[The determinant may also be found by another rule. Taking 
as before only the terms of the second order in the Lagrangian 
function L = T—V, let us separate the terms of the form 

Q = Cje" + CJ<I>' + CJ'ff> + &a 

In the remaining terms put ^ = ^m V— 1, ^' = ^m V— 1, and 
so on, and write down the discriminant If the system oscillates 
about a position of equilibrium the terms represented by Q are 
absent and the discriminant thus formed will be the determinantal 
equation giving the required values of m. But if the system 
oscillate about a state of steady motion we must modify the dis- 
criminant by adding some quantity derived from Q to each term. 
To find this, write above the columns ff, <f>, &c. and before the 
rows 0, if}, &c. Consider any term, say the term in the column 
^' and row 0. We must add to that term (Cjj— G^ m, where 
(7,g — Cj^ is the excess of the coefficient of ^0 above the coefficient 
of ^ff m the expression for Q, Since the determinant is un- 
changed by writing — m for m, we may, if preferred, add the 
excess of the coefficient of 0^ above the coefficient of 0^ provided 
we adhere to one order throughout.] 

6. If in the determinantal equation we write — m for m, the 
rows of the new determinant are the same as the columns of the 
old, so that the determinant is unaltered. When expanded, we 
shall have an equation which contains only even powers of m. 

The condition of dynamical stability is that the roots of this 
equation should be all of the form 

If we write m' = — Jp^ the roots of the transformed equation 
must be all real. 

In the case of equal roots of the form m= + /8V— 1, it has 
been shown in Art. 5 of Chap. i. that it is necessary for stability 
that the proper number of minors in this determinant should 
vanish. If there be two equal roots, these roots must make all 
the first minors zero ; if three equal roots, all the first and second 
minors must vanish, and so on. 

7. When the system depends on many co-ordinates the labour 
of expanding this determinant is often considerable. Methods of 
evading this in certain cases will be given in the following chap- 

R. A. 4 
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ters. But when the development is necessary we may proceed in 
the following manner. Let the determinant be written 



■SxtWl'+Ou 


B„m^ + a„ 


AX + «i. 




+ i^„m 


+ F„m 


B,X + o„ 


B„m' + a„ 


■S„m' + a„ 


+ F^m 




+ F„m 



&c. 



&c. 



&c. 



and let 



^=^110 + B^0<f)'\- ,.., 



r» .^ IT 

We- know that the determinant when expanded is of an even 
order, hence all odd powers of m must finally vanish. Let us 
expand the determinant in powers of the F's, The first term is 
the discriminant of Bm^ + a, this term is independent of the F^s. 
The terms which contain the first powers of the jF's are obtained 
by erasing any one line of this discriminant and replacing it by 
the corresponding F terms. But these terms all vanish and we 
need not describe them minutely. The terms containing the 
products and squares of the ^'s may be obtained by erasing every 
two lines of the discriminant and replacing them by the corre- 
sponding jP terms. Thus if we erase the two first lines we have 
the determinant 






F,,m 


F^m &c 


-F.^ 





F„m &c. 


B„m* + a„ 


B^m' + a„ 


B^m' + a^ &c. 


&c. 


&c. 


&c. 



and so on for all the other rows taken two and two. The terms 
which contain the cubes and all odd powers of the -Ps vanish, 
while the terms which contain the fourth powers may be obtained 
by erasing four lines of thh discriminant and replacing them by 
the corresponding jPs. 

When the determinant has been expanded, we have an equa- 
tion of an even order to find the values of m. We may therefore 
employ the short method of Art. 5, Chap, in., to obtain the 
Sturmian functions. 

8. Necessary and sufl&cient tests of the stability of the motion 
of a system of bodies are given in the preceding pages. But it 
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is assumed, as explained in Art. 2, that the co-ordinates have been 
properly chosen. They are supposed to have been so chosen that 
the coefficients in the expanded Lagrangian function are all 
constants. When this is not the case we must discover the 
proper co-ordinates to which the system must be referred before 
we can apply the test of stability. But when the motion is steady 
this is not difficult. 

There are obviously many such systems of co-ordinates, and 
one set may generally be found by a simple examination of the 
steady motion. If there are any quantities which are constant 
during the steady motion, they will often serve for some of the 
co-ordinates. Others may be found by considering what quan- 
tities appear only as differential coefficients or velocities. Practi- 
cally these will be the most convenient methods of discovering 
proper co-ordinates, since no further change will then be necessary 
and we may at once form the determinant of stability. But if 
these methods fail we may adopt the following analytical method 
of transforming (where possible) the general Lagrangian function 
with variable coefficients into one with constant coefficients. 

9. Let the Lagrangian function be 

+ ^A^^0' + Aj9<f>+... 
+ i B J* + B,,ff <!>' + ... 

where the coefficients are all functions of t and the co-ordinates 
0, if), &c. have been so chosen as to vanish along the steady motion. 
We have therefore for the steady motion 

The oscillations about the steady motion ate given by the 
terms of the second order. Our present object is to transform 
these to others with constant coefficients by the following substi- 
tutions : 

=p^x +PJ/ +p^z -f &c., 

&c. =&c., 
where the p's, q\ &c. are functions of the time at our disposal. 

Substituting and equating the coefficients of x^, y", &c. to 
unity, we have as many equations of the form 

^B,y + ^B„^ + B,,pq+.., = l (1) 

4—2 
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as there are co-ordinates. Equating the coefficients of the pro- 

n — 1 
ducts xy\ ojV, &c. to zero we get n — 5— equations of the form 

Ai;^ii>.+-SM?i?2+^i2(M«+Mi) + - = o (2), 

supposing that there are n co-ordinates. 

Equating to the constants a^, cTj, ... the coefficients of xx', yy\ 
&c. having subtracted the differential coefficients of (1) we have 
n equations of the form 

Adding the coefficients of xy' and x'y and subtracting the 

w — 1 - . 
differential coefficients of (2) we have n — ^ — equations of the 

form 

Equations (1), (2) and (4) give n* equations to find the n* 
quantities jPjjPj, &c., q^q^, &c. The solution of these equations is a 
purely geometrical problem. If we construct the two quadrics 

^B,,0' + ^B^<f>'-\-BJ<f>+.., = 1, 

and refer them to their common conjugate diameters, by writing 

0=PiX+p^-\-... 
^ = qiX + qj/ '\- ... 
&c. = &c., 

making the first quadric to become what we may call a sphere by 
projection; the values o{pj)^ &c., q^q y &c. thus found are the values 
required to make some of the coefficients in the Lagrangian 
function become constant. These values must of course make all 
the other coefficients of the second order in the Lagrangian function 
constants also, and thus we haven(n+]) analytical conditions 
that the motion should be steady. 

It might be supposed that greater generality would be obtained 
by replacing the zero's of equations (2) and (4) or the unities of 
(1) by arbitrary constants. This may be convenient in practice, 
but as we know that by a subsequent real change with constant 
values of p^p^, &c., q^q^, &c., we can render them zero or unity, it 
simplifies the argument to perform the two transformations at 
once. 
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10. The geometrical problem just alluded to admits of a real 
solution whenever one quadric can be projected by a real projection 
into a sphere. The problem then becomes that of finding the 
principal axes of the other. This is just our case, since the ex- 
pression 

i5,^» + /3,/f 4... 
is necessarily positive for all values of &^ , 

It is unnecessary to describe here the mode of solving this 
problem. It is sufficient to say that it may be reduced to the 
solution of the symmetrical determinantal equation 






= 0, 



where Du, D^, &c. are the coefficients of 0^, 5^, &c. in the second 
of the quadrics. The roots of this equation are known to be real 
when the suppositions just mentioned are satisfied. 

11. In order to examine the fundamental determinant in 
Art 4 a little more closely, let us suppose it reduced to depend on 
three co-ordinates. We may then have the advantage of a geo- 
metrical analogy. Let the co-ordinates be f r), f, and let the 
equations of motion be written 

^'"^-'■H -oi H .4 n 



+»a ' ^ --^s 




j,+(a.|-^„){=o. 



Let a geometrical point P move in space so that its co-or- 
dinates referred to any axes are f, 17, f. Then the position and 
motion of the point exactly give us the position and motion of 
the system. 

12. Looking at the equations of motion just written down we 
see that they are similar to those which give the oscillations of a 
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system about "a position of equilibrium, but that there are in 

addition terms E ^, F-j, &c. Th^ general effect of these terms, 

at at 

as will appear from what follows in the subsequent .chapter, is to 

increase the stability. If we transpose these terms to the other 

sides of the equations, we may regard them as impressed forces 

acting on the system, whose resolved parts in the directions of the 

axes f , 7j, f, are 

at at 



We see at once that 



^di'^-^dt 



HX+FY+OZ- 



dt dt dt 








I < 



so that these forces are at once orthogonal to the path of the 
representative point P and also orthogonal to the straight line 
whose direction cosines are proportional to JE, F, 0, These forces 
are therefore of the nature of centrifugal forces^ as if they were 
produced by the rotation of the system about this straight line. 

13. We may show that the straight line (J?, F, Q) is fixed in 
space. To prove this, let us transform our co-ordinates from f,i7, f 
to a?, y, fir, where a?, y, z are connected with f , 17, J by any linear 
relations, such as 

^^a^x^-hjf + c^z] 

t^^a^x-k-ljj + c^zl 

Let the portion of the Lagrangian function under consideration 
(Art. 1) be 



then 



^= ^12"" ^Sl> -^= ^28 " ^82> -^^ ^81 "" ^13- 



Substituting for f , 1;, ? their values in terms of a?, y, z, we find that 
the difference between the coefficients of xy' and xy is 



«.«, 


+ E 


«. a» 


+ F 


Oil«i 


\\K 




KK 




hh 
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with similar equations for E' and F\ But if /* be the determi- 
nant of transformation 



«1«. 


+ f 


«»«« 


+ 1? 


o.«. 


h,\ 




\\ 




J, Ji 



with similar equations for x and y. The ra<zo5 of ^, -P^ G^ are 
therefore transformed as if they were co-ordinates. If the trans- 
formation be a real transformation of Cartesian co-ordinates, let 
lengths each equal to unity be measured from the origin along the 
axes Of, O17, 0£ thus forming a tetrahedron whose volume is F. 
Let a similar construction be made for the new axes, forming a 

F 
tetrahedron of volume F. Then* ^^-y.. Hence the quantities 

=, =, =^ may be transformed as if they were lengths measured 

-pt -pt Qf 

along the axes and become ^^ , ^77 > T77 • If both systems of co- 

1 1 

ordinates are rectangular we have F= « , ^ = 5 • 

E F Q- 

Let (o be the resultant of j^, y.^ p. , then to may he regarded 

as a fixed length measured from the origin along a straight line 
fixed in space. 

Let V be the velocity of the representative particle, the angle 
between the direction of this velocity and the axis whose direction 
cosines are proportional to E, F, O. Then the resultant of the 
forces X, F, Z is easily seen to be 2v Fo) sin acting perpendicular 
to the axis and to the direction of the motion. We might call the 
straight line (EFG) the axis of the centrifugal forces. 

[* Let {^iVil^i)> tta^yafa)? (^8% ft) ^^ *^® co-ordinates of three points -4, jB, (7 re- 
ferred to any oblique co-ordinates. Let us find the volnme V of the tetrahedron 
of which these and the origin are the angular points. Since the volume vanishes 
when any angular point as C Hes in the plane containing the origin and the other 
two A, Bf the expression for the volume must contain the factor 



A*= 



^1 ^n ^3 

""Ji 'fi 'fz 
f 1 fa f 8 



The volume is evidently an integral rational function of the co-ordinates when 
the axes are rectangular and the plane AOB is taken as the plane of xy, it 
easily follows that this is true for all axes. Since this function cannot he of more 
than the third order, we have V'=Mfi, where M is independent of the co-ordinates 
of At Bi C, When the points A^ B^ C are on the axes at unit distances from the 
origin, let Fhe the volume of the tetrahedron. In this case ai=1, and .*. M= V. 
We have therefore in all cases F' = Vfi, 

In the text, let the extremities of the unit lengths measured along the axes of 
as, y, z he called A, B, C, Then the (^vi) co-ordinates of A, B, C are (a,, a«, Og), 
(&i» &a» ^8)» (^i> ^v <^s)> respectively. Hence hy what has just heen said k'=5K>.] 
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14. The expressions for the co-ordinates in terms of the time 
will in general contain as many periodic functions as there are 
co-ordinates. If the initial conditions are such that each contains 
one and the same periodic function, the motion recurs after a con- 
stant interval and the system is said to be performing a simple or 
harmonic oscillation. 

If the system be oscillating about a position of equilibrium, 
with a Lagrangian function 

we know* that the harmonic oscillations are represented hy rectili- 
near motions of the representative particle, and that these are 
along the common conjugate diameters of the two quadrics ^ ' 
equations are 

•»n 1 + ^12^7 + ...=6 

where a and b are two constants chosen to make the quadrics real. 
Let us consider what are the harmonic paths of the representative 
point when the system is oscillating about a state of steady motion. 

In any harmonic vibration we have f = i cos (Xi + a) * with 
similar equations for rj and f. Hence 

df^^^^' de — ^' 5? — ^^• 

Substitute these in the equations of Art. 11. DiflFerentiate and 
substitute again. Multiply by f, r), f, add and integrate, we obtain 

where c is some constant. The harmonic path lies on this quadric, 
which has a common set of conjugate diameters with the two 
quadrics a and b. 

If we resume the result of the substitution of -p , &c. in the 

equations of Art. 11, and multiply by -&, F, respectively and add, 
we obtain 

[(5^^+^i,F + 5^(?)| + i&c.]V + [(^^^+...)? + &c.]=0, 

* A short paragraph in Thomson and Tait*s Natural Philosophy ^ page 273, is 
the only notice of this which the author has discovered. 
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which is a plane, and is diametral to the straight line {EFO) with 
regard to the quadric c. 

The harmonic paths are therefore ellipses. The three harmonic 
planes are diametral to the same straight line and this straight 
line is fixed in space, being the axis of the centrifugal forces. 

If we eliminate \ between the equations to the plane and the 
quadric c, we get a cubic surface on which the three harmonic 
conies lie. 

If Ey jF, are zero, which is the case when the system oscillates 
about a position of equilibrium, the quadric c becomes a cylinder. 
This mafbe conveniently showA by referring the system to such 
co-ordinates that the coefficients B^, By^, B^, A^y A^, A„ are all 
zero. In this case the diametral plane of every straight line passes 
through the axis of the cylinder. The harmonic oscillations are 
therefore rectilinear. 

If B be the length of that semidiameter of the quadric (c) 
which is parallel to the fixed straight line (-iS', F, G), it may be 
shown that the 

Product of the axes ) R 2c 

of the quadric c J JW^TF^^O^ ' ^ * 

If E, F, G are all zero, and their ratio is indeterminate, B is 
any diameter. Hence one of the axes of the quadric (c) must be 
infinite and the quadric will be a cylinder. 

[If the quadric (c) be a cylinder and E, F, O are not all zero, 
we must have either \ zero or B infinite. In the latter case the 
axis of the cylinder will coincide with the axis of the centrifugal 
forces.] 

The quadric (c) has also the following geometrical property. 
Let the lengths of semidiameters of the quadrics (a) and (6) drawn 
parallel to the axis of the centrifugal forces be p and p. Through 
the intersection of these quadrics describe a quadric so that the 



Product of) 2\fab 

its axes 

- ' " ' ^ A/ I 

This quadric is similar to the quadric (c). 



iv'-?)^ 






15. The introduction of the representative point to exhibit 
the motion of a system may appear somewhat artificial. If how- 
ever we properly choose the co-ordinates the particle moves exactly 
as a free particle, and we might reduce the problem of finding the 
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oscillations of a system to a problem in Dynamics of a particle. 
Refer the quadric 

to its principal axes and let the equation thus changed be 

B^ ^ + -Sffl' -- + ... = 6. 

Since J?u , B^^ &c. are positive quantities, we may put 

V^'?i=a?, *jBj7)^ = y, &c. 

The quadric has thus been "projected" into a sphere. Let a?, y, z 
be now chosen as the co-ordinates of the system and let the 
Lagrangian function be expressed in the form 

the terms of the first degree being omitted as not necessary to our 
present purpose. The three equations of motion at the beginning 
of Art. 11 take the form 

&c. = 0, 

which are the three equations of motion of a free single particle 
whose co-ordinates are x, y, z under the action of forces whose force 
function is 

and a force acting perpendicular to the path and also perpendicu- 
lar to a fixed straight Une, the force being proportional to the 
velocity. 

16. As an illustration of this theory, let us here make a short 
digression. However the particles of light may oscillate, whether 
in a rotatory or linear manner, we know the motion is related to a 
certain plane called the plane of polarization. It may be shown 
that any harmonic oscillation about a position of eqwilihrium may 
be represented by a rectilinear oscillation of the representative 
particle. Let us represent the motion at any point of the ether 
by a rectilinear oscillation in a direction perpendicular to the 
plane of polarization. This would be FresneFs Vibration. The 
representative particle, as just shown, would not necessarily move 
as if it were a free single particle. But let us assume (and a proof 
is not necessary to bur present purpose) that when the oscillation 
is drawn as above described the motion in the plane of the front 
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is the same as that of a free particle, while that perpendicular is 
not free. On this assumption we see that Fresnel m his theory 
of double refraction is justified in taking actual instead of relative 
displacements, for it is the representative particle he is considering. 
He also neglects the force normal to the front, for the particle 
moves as a free particle only in the plane of the front. These 
general remarks are not meant to explain Fresnel's theory, but 
merely to show how the representative particle may be used to 
replace a complicated motion. 

17. [When a system is performing a harmonic oscillation 
about a state of steady motion or about a position of equilibrium, 
the motion repeats itself continually at a constant period, that is 
to say, the values of the co-ordinates recur at this interval. This 
is the chief peculiarity of a harmonic oscillation. When the 
oscillation is about a position of equilibrium, the representative 
particle oscillates in a straight line whose middle point repre- 
sents the position of equilibrium. Thus the system passes through 
the position of equilibrium twice in each complete oscillation. 
When the oscillation is about a state of steady motion the path 
of the representative particle is an ellipse whose centre is at the 
point occupied by the system in steady motion at the same in- 
stant. Thus the system does not in general ever coincide with 
the simultaneous position of the system in the undisturbed or 
steady motion. When a system is disturbed by a small impulse 
fropi a state of steady motion, it will in general describe a com- 
pound oscillation made up of at least two harmonic oscillations, 
at the instant of disturbance these two neutralize each other so 
that in the disturbed and undisturbed motions two simultaneous 
positions are coincident. But it is clear that this cannot occur 
again unless either the periods of the two harmonics are com- 
mensurable or the period of one of them is infinite.] 

18. [In some cases the ellipse degenerates into a straight line. 
Thus if the quadric (c) be a cylinder the diametral plane of the 
axis of the centrifugal forces will pass through the axis of the 
cylinder, and thus the harmonic oscillation corres;ponding to this 
particular value of \ will be rectilinear. In this case the system 
twice in each oscillation passes through the position it would have 
occupied at the same instant in the undisturbed motion. 

The quadric (c) has a common set of conjugate diameters with 
the quadncs (a) and (J). Hence if (c) be a cylinder, its axis must 
be parallel to one of the three common conjugate diameters of (a) 
ana (J). If we refer the quadrics (a) and (b) to their common 
conjugate diameters, they take the form 

^n'r+5«v+-B„r=26r 



60 THE MODIFIED [CHAP. 

The cylinder which passes through their intersection and has 
its axis parallel to the diameter fis found by eliminating f * between 
these equations. We see therefore that BJ\^ + A J = 0. If then 
the axis of the cylinder cut the q^uadrics (a) and (6) in D and 1/ 
respectively, we find that for this oscillation 



, oiy' -a 



^ ~ OD^' b ' 

It has already been shown that when this value is finite, the 
direction of OLD' is along the axis of the centrifugal forces.] 

19. [In some cases two or more of the values of \ are zero. 
In these cases the co-ordinates will have terms of the form nt + €, 
where n and e are two small constants. When, as explained in 
Art. 3 of this Chapter, there are several parallel states of steady 
motion, these terms imply that the motion is stable about a state 
of steady motion very nearly the same as the undisturbed motion 
but not coincident with it. The actual undisturbed motion, unless 
n is zero, is unstable in the sense that if a proper disturbance 
be given to the system, the system will depart widely from the 
positions it would have simultaneously occupied in the undisturbed 
motion.] 

20. In many cases of small oscillations it will be found that 
the Lagrangian function T—V is not a function of some of the 
co-ordinates as 0, <f>, &c. though it is a function of their differential 
coeflScients 0', <^', &c. In such cases the steady motion will be 
usually given by constant values of these differential coeflScients, 
while the other co-ordinates as f, rj, &c. are also constant. It is 
evident that the determinantal equation of Art. 4 is needlessly 
complicated. It is clear that there will be as many pairs of roots 
equal to zero as there are co-ordinates 0, <j>, &c. It will be an 
advantage to eliminate 0", <f>', &c. altogether from the Lagrangian 
function, and to find the remaining roots by operating only with 
the co-ordinates ^, rj, &c. We shall thus obtain a determinant with 
just as many rows as there are co-ordinates of the kind f , rj, &c. 

Let i/j be the Lagrangian function expressed as a function of 
ff<l>, &c. ^, ^'^', &c. Let L^ be its value when 0'(l> are elimi- 
nated, so that i/g is a function of ^, ^rj', &c. only. To effect this 
elimination we have the integrals 

dT_ dT_ 

where c^yC^y &c. are constants. Then 

dL,_dL dL, dff dL, d£ 
rff ~ d^ "^ dff "rff '^W C 
_di, dff #' 

~ df '•'''' 3F "*■"* df 



IV.] LAGEANGIAN FUNCTION. 61 

dZ/j _ dL^ do d^ 



d^dI^_dL,_d^dL^_dL, . (dd& d0' 
dt df df ~ di df df 



u UJU^ axyj _ t* UJJ^ »xy, / d d^ dff\ a 



But ^-.^-%} = 0. 



dt d^' d^ 



Hence if we take 



U = L-c/-c,<^'-cfcc. 

and eliminate ff, <\> by help of the integrals -[^ = Cj, ^ w = Cj, 

we may treat L' just as we do the Lagrangian functicm L. The 
equations giving the small oscillations about the steady miction 
vnll be. 

d dli dL' f\ J A 

dtdr""d|"=^' '^^•"^• 

The function L^ may be called the modified Lagrangian fv/nction. 

[It should be noticed that this is equivalent to a partial use of 
Hamilton's transformation of Lagrange's equations. Sir W. R. 
Hamilton eliminates all the differential coeflScients 0", 0', &c. by 

jrp jj rp 

the help of equations of the form ^ = w, ^ = v, &c. where u, v, 
&c. are made to be new variables*. In our transformation only 

* The Hamiltonian transformation of Lagrange's equations bears a remarkable 
analogy to the transformation of Beciprocation in Geometry. This may be shown 
in the following manner. 

When the system has three co-ordinates 0, <f>, ^, we may regard 0^, </>\ y// as the 
Cartesian co-ordinates of a representative point P. The position and path of P will 
exhibit to the eye and will determine the motion of the system. Let u, v, ti; be the 
Hamiltonian variables, so that 

dT. dt. dT. 

where T, is the semi vis viva expressed as a function of 0, <f>, ^, ff, <f>', yj/. Then 
u, Vy w may be regarded as the co-ordinates of another point Q whose position and 
path will also determine the motion of the system. 

If the semi vis viva be given by the general expression 

Ti = J Jii^+-4i2dV+ ... 

it is dear that the point P always lies on the quadric T^=U where Ui^ the force 
function and the co-ordinates 0^ </>, ^ have their instantaneous values. The point 
Q must therefore lie on another quadric which is the polar reciprocal of the first 
with regard to a sphere whose centre is at the origin and whose radius is equal to 
1J2U, The equation to the reciprocal quadric is therefore 

u V w =Ut 



^2- 2A 



U Ay^i Jjj ^j3 

w Aj^ A^^ A^ 
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those new variables are introduced which would be constants in 
Sir W. R. Hamilton's transformation. 

This remark suggests an extension of the process. If i be a 
function of ^, <f>, &c. as well as of ^, d>, &c. the quantities c^, c^, &c. 
will not be constants. We express tnis by writing u, v, &c. instead 
of Cj, c,, &c. Suppose we wish to eliminate some of the differential 
coefficients, viz. ff, (f), &c. and to retain the remaining ones, viz. 
^, 7/, &c. If we put 

i' = i-w5'-t;<^'-&c. 

we may easily show as in the preceding page that 

d dL dJj /v o /v 

where A is the determinant, called the discriminant, which may be formed from 
the determinant just written down by omitting the first row and the first eolnmn. 

This is a general expression for the Hamiltonian function and agrees with that 
which may be deduced from the result in Art. 21, when all the yariables are trans- 
formed by the Hamiltonian process. 

Since the polar reciprocal of thtf polar reciprocal is the original quadric, it 
follows that 

ffJ''^^ ^'^^ ^^^T^ 

^'dH' ^-d;"' ^-d^' 

which are three of the six Hamiltonian equations. 

We may also show geometricaUy that if the coefficients of T^ be functions of 

dT dT 
any quantity ^, then -^= - -^^ • To prove this we notice that if as, y, x; be the 

co-ordinates of a point P, situated on a radius vector OP' of a quadrio 

jyp,) •• The quadrics 

Ti = l and T^=l may be regarded as polar reciprocals of each other with regard 

to a sphere whose radius is 1J2 and whose centre is the common centre of the 
two quadrics. Let P be any point on the quadric 7^=1, and let the radius vector 
be produced to Z so that OP . 0Z=% then the quadric 2",=! touches a plane 
drawn through Z perpendicular to OP and Q is the point of contact. Let these 
quadrics be slightly altered in consequence of a variation of 0^ so that their 
equations are now Ti+dTj^=l and T^+dT^=l. Let OP and OQ produced cut 
these quadrics respectively m P^ and Q\ Then 



2'x+<in=(g>y. r.+<ir.=(^)'. 



Now if Z' be a point on OP produced so that OP' . OZ'=OP . OZ.ihe quadric 
T^+dT^mH touch the plane drawn through Z^ in some point q near Q'. The 
point Q' will therefore lie very nearly in the tangent plane, so that by similar 
triangles 

OQOZ _ OP' 

OQ: " OZ' '" OP ' 

Since each of these ratios is indefinitely nearly equal to unity, it follows that 
dTi=-dT^, 

If we put L=Ti-{-U and H=Ta- V, Lagrange's equations may be written in 

the forms m'= -j^ , «'= jt» «^'=-rT« Hence we have 
dQ d(/> df 

, dH , dH , dK 

-~=5?» -^=d0' -'^^S^' 

which are the remaining three of the Hamiltonian equations. 
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We have thus as many equations of the Lagrangian form as 
there are variables f , 17, &c. Also since w = t^? , &c. we have by 
differentiation 

with similar equations for j>\ &c. By Lagrange's equations we 

obtain 

dL' , p 
^=w, &c. 

Thus we have as many sets of equations of the Hamiltonian 
form as there are variables ^, ^, &c.] 

21. We may eflfect this elimination once for all and find a 
definite expression for L\ 

Let the kinetic energy be 

Then the integrals used will be 

2flfl^ + 2fl^^' + ... = Ci— Tb^^ — T^i{ ^ ... 

T^ff '\-T^^ '\' ... = Cg— if^f — r^,i7 — ... 

&c. = &c. 

For the sake of brevity let us call the right-hand members of 
these equations Ci — -Z, Cj — F, &c. Since T is a homogeneous 
function, we have 

-i^(Ci-^-if(c,-r)-&c.J 

If we substitute in the second line the values of ff, ^', &c. found 
by solving the integrals just written down, we have 



i' = TfX'+2V,fi?' + &c.-F 



''"2A 



c,-Z c^-Y 
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where A is the discriminant of the terms in T, which contain only 
ff, <\>\ &c., and may be derived from the determinant just written 
down by omitting the first row and the first column. 

We may expand this determinant and write it in the form 



^'=^«Y+^f,rv+&c.- V 



+ 



2A 













t^ 



T, 



** 



•• ••• ••« 



+ 



2A 



X Y 

Y 1$^ 1^ 



1 

A 







X 

Tea 
Tm, 



T. 



^ 



• • • • • 



• • ^ • • • 



where X, T, &c, stand for 

&c. = &c. 

The first of these three determinants will contain only the 
constants Cj, Cg, &c., and the co-ordinates f , rj, &c. The second will 
not contain Ci, Cj, &c. but will be a quadratic function of f', 7j\ &c. 
The last determinant will contain terms of the form f, t) with 
variable coeflScients which may also be functions of c^, c,... 

Since ^^, 1717', &c. are all small quantities, it is clear that this 
expression for L' when expanded will take a form precisely similar 
to that given in Art. 2, only that we have fewer variables to deal 
with. 

22. As an example, let us consider the following problem. 

A body has a point O which is in one of the principal axes at 
the centre of gravity G fioced in space. The body is in steady 
motion rotating with angular velocity n about OG which is vertical. 
Find the conditions that the motion may be stable. 

Let OA, OB, 0(7 be the principal axes at and let OG co-. 
incide with the vertical OZ in steady motion. Let f , 17 be the 
direction cosines of the vertical OZ referred to OA, OB. Let 
«i, o)j, o)g be the angular velocities about the principal axes at 0. 
Then to the first order 

c»i = V + oyJ^X ^ 
©2 = - 1' + 0)3771 • 
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Let 6 be the angle ZOC, ^ the angle the plane ZO (7 makes 
with a plane ZOX fixed in space and ^ the angle it makes with 
the plane J. (9(7 fixed in the body. Then 

f = — sin^-cos^] 
rj= sin ^ sin <^J * 

t»3 = ^' + '^'<50S<9 

='^'+(x'-f){i-f)> 

putting ;^ = ^ ^ '^. \ye easily find 

If then A, £, Che the principal moments of inertia at 0, the 
Lagrangian function is 

+^('?'+xf)'+|(-r+x9)' 

where JIf is the mass of the body, «iid h = 00-, 

Since x ^ absent from the equation we have the integral 

dL_ 

which gives 

j^'^n + terms g( second order. 

Hence 

-const + ^," + |r 

+ {( J - 0) «• + Jf^A} 1+ {(5 - (?)«• + %M I* 

Using this as the Lagrangian function we easily find 

(A - 0)71* + Mgh - Bm', (A + B-C)nm ' =0. 

-(A+B-C) nm, {B -C)n* + Mgh - Am 

B. A. 



66 THE MODIFIED [CHAP. 

The roots of this equation to find m must for stability be of 

the form +/8V — 1. Putting m* = — \* we have a quadratic to 
find \'. The roots of this quadratic must be real and positive. 

If -4 = J?, as in the case of a top spinning with its axis vertical, 
we have 




- 4^AMgh\ 



2A 

The motion is stable or unstable according as C*n^ is greater or 
less than ^AMgh, If CV = 4iAMghy the equation has equal roots 
and as the first minors are not zero the motion is unstable. 

23. [As another example of the use of the modified Lagrangian 
function, let us consider a case discussed by Prof. Ball in the 
Notices of the Royal Astronomical Society for March, 1877. In a 
problem in Physical Astronomy, we want the relative co-ordinates 
of the system, while its absolute motion in space does not concern 
us. Lagrange's equations involve both the relative and absolute 
co-ordinates, and are therefore not particularly well adapted for 
such problems. By using the modified Lagrangian function, we 
may eliminate the absolute co-ordinates. 

Let the system have n co-ordinates, let us choose as three of 
them the co-ordinates of the centre of gravity of the whole system, 
viz. 0, (f), '^. There will then remain n — 3 co-ordinates which are 
independent of these.. Let T be the kinetic energy of the system 
relative to its centre of gravity, V the potential energy, M the 
whole mass, Their the Lagrangian function is 

L=^M{ff'+ ^'»+^'»)+ r- V. 

In problems in Physical Astronomy the potential energy is a 
function only of the relative positions of the bodies, and is there- 
fore independent of 0, <f), yft and their differential coefficients. We 
have therefore 

dL dL ^ dL _ 

dff~^^' df"^^ df*"^'' 

Hence the modified Lagrangian function is 

L' — T -^ F— a constant. 

It is independently clear that we might take this as the 
Lagrangian function, for the first three terms of L do not enter into 
any one of the Lagrangian equations, except the three formed by 
differentiating with regard to ff, <f>, '^\ 

The function T is made up of two parts, (1) the kinetic energies 
of the rotations of the bodies about their centres of gravity, which 
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we may call 7/, and (2) the relative kinetic energies of the several 
bodies, each coUected at its centre of gravity, which we may call 
T'. Let m., w , &c. be these masses; a?-, x^ &c. the abscissae of 
their centres of gravity referred to the centre of gravity of the 
whole as origin. Then, accents denoting diflFerential coefficients 
with regard to the time, we have 

m^xl -f m,aj,' + &c. = 0. 
Let us square this and write 

2xX = ^'" + xi^ - « - <)•. 

If we examine the coefficient of any power as a^'* we see that it 
is 

Hence the square becomes 

2w %7nx'^ — %m^m^{x( — a?/)* = 0. 

Similar expressions hold for the y and z co-ordinates. Hence 
on the whole we see that the relative kinetic energies of the several 
bodies collected at their respective centres of gravity is 



22m ' 

where v is the relative velocity of the centres of gravity of the 
masses m^, m^. If we express this in any kind of co-ordinates, we 
may use the Lagrangian function L' to find the relative motion. 

The expression for 7/ agrees with that given by Prof. Ball, but 
his demonstration is quite diflFerent. Prof. Cayley has given 
another demonstration in the same number of the Astronomical 
Notices. 

The Lagrangian function thus found may be still further 
"modified." To avoid symbols of summation, let us consider the 
case of three particles moving in one plane under their mutual 
attractions. Let the separate masses be m^, m^ m^, and let /i be 
their sum. Referring the system to m^ as a central mass, let the 
distances of tw,, m, from m^ be respectively r ajid p, and let the 
opposite side of the triangle be It, Let the interior angle be- 
tween r and p be <^ and the exterior angle between r and R be 
X' Let 6 be the angle r makes with some fixed straight line in 
space. We easily find 

where A = m^m^ + m^j? + m^m^I^y 

C = zr- Xmrnv^y 



^1 



¥ 
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and V is the relative velocity of the masses m, m' calculated on 
the supposition that mi is fixed, and that the straight line r has 
no rotation round wii. Thus A, JB, C are all functions of r, p, <f> 
and their diflferential coefficients with regard to the time. 

If we only want the changes in the form and magnitude of 
the triangle joining the three particles, we may eliminate 0" by 
means of the equation 

Aff+B = c. 

We then find as our modified Lagrangian function 
which contains only the three co-ordinates r, p and ^]. 

« 

24. When the geometrical equations contain differential 
coefficients of the co-ordinates ^, ^, f, &c. of the system with regard 
to the time, we cannot express the co-ordinates a;, y, z of any 
element of a body in terms of f , rj, f, &c. by means of equations of 
the form 

y =/2(?> ^> £ &c., t\ 

It follows, as is pointed out in our books on Rigid Dynamics, 
that Lagrange's equations cannot be employed in the form 

d dT_dT__dV 
dt d^ df df ' 

In many of the most interesting problems in Rigid Dynamics, 

di 
it so happens that the geometrical equations do contain -~, 

dn 

-— , &c. For example, let a sphere be set rotating about a vertical 

diameter and be on the summit of a perfectly rough surface of 
any form. If a small disturbance be now given to it, the sphere 
may roll round and round the summit. During this motion the 
velocity of the point of contact is zero, and our mode of re- 
presenting this analytically in terms of the co-ordinates will give us 
two equations of the form 

To include such cases the equations of motion must be modi- 
fied. If L be the difference between the kinetic and potential 
energies, all the Lagrangian equations may be written in the form 

fd dL dL\ ^^ . fddL dL\ ^ , « ^ 
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where Sf , Srj, &c. are any small arbitrary displacements consistent 
with the geometrical equations. But if these geometrical equa- 
tions be given in the form 

&c. = 0- 



these arbitrary displacements must satisfy 

&c 



. = 0) 

.=or 



K this were not the case, the geometrical displacement of the 
body given in applying Virtual Velocities would not be such as 'to 
cause the unknown frictional forces, &c. to disappear. Using the 
principle of Indeterminate Multipliers, we get 

d dL dL ^ j^ TT r\ 



&c. = 






These joined to the geometrical equations 

5.=or 



&c. 



are suflScient to determine the unknown co-ordinates f, 17, &c. and 
the multipliers \ fi, &c. 

It will be more convenient to write these equations in the form 
±dL_dL dG dH 

ddL dL dG dH 

&c. = 0, 
the geometrical equations being 

(? = 0, 5^=0, &c. 

It is of course obvious that these indeterminate coefficients 
X, fji,, &c. are merely the frictions or other resistances introduced 
into the equations in a convenient form. 



25. In order to apply these equations to the oscillations of a 
system about a state of steady motion, it will be convenient to 
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change the co-ordinates f , % &c. into others ^, ^, &c. which vanish in 
the steady motion. Let L be thus expanded in powers of 0, <^, &c. 
as explained in Art. 1, and let P and Q have the meaning given 
to them in Art. 5. 

Let us then put 

^ = a + ^, V = i9 + f, &c. 

\=\-f\, /i = A*o + A*i> &c. 

where a, /8, \, /Lt^,, &c. are the. values of ^, rf, \, /it, &c. in steady 
motion. The geometrical equations will* then take the form 

(?=<?, (a+^ + G,(/8+f)+&c. = 0| 

&c. = 

and the equations connecting S^, S^, &c. will be 



(?,S^ + (y,S0 + &c. = 

&c 



.=0) 

.=or 



In these equations O^, 0^, &c. are functions of 0, ^, &c. Let 
a square bracket indicate that the value of the inscribed quantity 
in steady motion is to be taken. Thus [GJ means the value of G^ 
when 0, <\>y &c. have all been put zero. 

The equations of steady motion may then, exactly as in Art. 3, 
be written 

&c. = OJ 

From these the relations which exist between a, ^, &c., \y fif^y &c. 
may be found. 

The equations of the oscillatory motion may be written 

^Id ^W ^"^^ ^'\dff\ ^'^^ '''ld&\ '^^ "' 

with similar equations for ^, '^, &c. 

« 

26. The final determinant written for two variables ^, ^, and 
two geometrical equations O and H in the notation of Art. 19, 
will be 



► 
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B,y - A,, + E,, 

- -do; , 'dH; 
'^'^»\_ddy'*'[d0] 


B„m*-A^ + E„ 


[dCn [dHl 
Idffj Idd'j 

'dcr 'du 

\_d^ \ \_d<\>' J 





B,,m* - A,, + E,, 




VdO'\ . rdGH 

[de\^Ydd'r 




-dHT , 'dH- 





-0. 



It will be noticed how very much this determinant is simplified 
if the values of \, /* in steady motion are zero. 

27. Let us apply these equations to the solution of the follow- 
ing problem. 

A heavy sphere rotating about a vertical diameter rests in equi- 
lihrium on the summit of a perfectly rough surface and being 
slightly disturbed makes small oscillations, find the periods. 

As the sphere moves about, its centre always lies on a surface 
which may be called parallel to the given surface. Let the high- 
est point of this surface be taken as the origin and let the axes of 
OS and y be the tangents to its lines of curvature at 0, so that the 
equation to the surface in the neighbourhood of is 



\pi PJ 



Let P be the centre of the sphere, PC that diameter which is 
vertical when the sphere is in equilibrium on the summit. Let 
PA, PB be two other diameters forming with PC a system of 
rectangular axes fixed in the sphere. Let the inclination of PC 
to the axis of Z, which is vertical, be 0, and let the vertical plane 
through PC make with the plane xz an angle '^y and with the 
plane GPA an angle ^. The vis viva 2T of the sphere will then 
be given by 

T^ i (^" + y") + i *' {(f + 1' COS ey + &' + sin^^.^'^}. 
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Let sin 5 COS -4^ = ^, sin 5 sin -^^ = 17, tten we have to the 
necessary degree of approximation 



^» + sin*5^'' = f'» + i7'M' 



Also let 6 H- "4^ = ;^. These transformations of co-ordinates are 
all permissible, because they do not involve any differential coeffi- 
cients with regard to the time. We thus find if i be the differ- 
ence between the kinetic and potential energies 

If a>g, ©y, CO, are the angular velocities of the sphere about 
parallels to the axes, the geometrical conditions are 



a?'" — a [©y,— o), ~ ] = 



where a is the radius of the sphere. These equations by well- 
known rules reduce to 



X 



Ir^'sin'^sin^ + ^cos.'^lr — ('^' H-^'cos^} — = 



a 






- ^ - ^' cos ^ sin ^ + fi' sin ^ + (^' + ^' cos 0)^ = 
& Pi 

expressing these in terms of our new co-ordinatea we have 



a 



rl 



a 



The position of the system has now been expressed in terms 
of such co-ordinates, that the coefficients in the governing func- 
tions L, G, H are all constant. See Art. 2. 

The steady motion is given by x, y, f , 17 all zero, and ^ = n. 
To find \^ fiQ we may use the equations of steady motion 

where q stands for any one of the co-ordinates. Taking j =s a? and 
J = y, we see that X^ = and fjL^ = 0. 

To find the oscillation we may substitute in the determinant 
and thus form the equation which gives the periods. As most of 
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the constituents of the determinants are zero, it will be more con- 
venient to form each equation directly from the standard formula 

where q stands for any one of the co-ordinates. Taking q in turn 
to be X, y, -y^ f , t) we find 



y -9 



a 



• 1 



;6»x" = 

A:* (V'-XV) + /*. = » 

Putting x = w these with the two geometrical equations are all 
linear and ready for elimination. 

Eliminating \, fi^ we have 



r-HnV-fSa."-f^^^=0 



Jfc* 






' 



Substituting for f , 17 from the two geometrical equations, we 
have 

^Pi <^P2 



a 



J(^ + a* „^ y _ . wi* / 



a 



— y -9 



a'p. 



X 



* • 



To solve these put x = Xcos {pt + j) y = Fsin {pt + q) so that 
p is the quantity required. We obviously have 

which is a quadratic to find p^. 

If pp p, have opposite signs the roots cannot be real, and the 
steady motion, must be unstable. If pj, p^ are both positive, so 
that the sphere is on the summit^ the motion is stable only if 



n*> 



CHAPTER V. 



Certain suhsidia/ri/ determinants a/re formsd from the dynamical de- 
terminanty <md it is shown that there mvst be at least as many 
roots indicating stability a» there a/re variations of sign lost in 
these subsidiary determinants, amd must eocceed the nwmher lost by 
am even number. Arts, 1 — 5, and 9. 

l^his is equivalent to a maximumi and minimwm criterion of stability 
with similar limitations. Arts, 6 — 8. 

Effect of eqvKil roots on this test of the stability of the system. Arts, 
10, 11. 

Examfiple, Art. 12, 

1, In order to test the nature of the roots of the determi- 
nantal equation, let us apply a method analogous to that by which 
Dr Salmon in his Higher Algebra proves the reality of the roots 
of the equation which occurs in the determination of the secular 
inequalities of the planets. 

2. Let A be the determinant which forms the left-hand side 
of the fundamental equation, let A^ be the determinant formed 
by omitting the p^ row and <f^ column. Let A^ be the determi- 
nant formed by omitting the first r rows and r columns. Thus 
Aj = Ajj. We then have by a known theorem in determinants 

AA =A A -^A A . 

It has already been noticed that if we change m into — m, the 
determinant A is changed into another determinant whose columns 
and rows are the same as the rows and columns of the first de- 
terminant. It easily follows that the minor — A^^ is changed into 
the minor -- A^^ by changing m into — m. Hence if 



Aj2 = <t> (m*) + w-^ (m*) 
then ^n — ^ (^^ — ^'^ (^') 



}• 
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Hence the product Aj^Aj^ is necessarily positive for ail negative 
values of m*. 

It also follows that if A-, vanishes for any negative value of w* 
then Ajj also vanishes for the ^me value of m*. 

3. When the determinant A^ vanishes, we have 

AA, = - A^„, 

so that A and A, must have opposite signs, or one of them must 
be zero. Consider then the series of determinants 

A, A,, A„ A3, 

each one being formed from the preceding by erasing the first 
row and the first column. We thus have a series of functions 
of m* whose degrees regularly diminish from the w*** to the first. 
As we may suppose the determinant A to have a row and a column 
of zeros added on at the bottom and right-hand side, but with any 
positive constant in the right-hand bottom comer, we may add 
to this series of determinants any positive constant. We have just 
proved that if any determinant of this series vanish for a nega- 
tive value of m', the two determinants on each side have opposite 
signs. The case in which two successive determinants vanish for 
the same value of m* will be considered afterwards. 

We may then use these determinants lA a manner somewhat 
similar to that in which we use Sturm's functions, provided no two 
successive functions vanish for the same negative value of m*. 
No variation of sign can be lost as we pass from m' = — oo to 
m* = except by the vanishing of the determinant A at the head 
of the series. And when a variation of sign is lost, it will be 
regained again at the next root, unless a root of the determinant A^ 
separates the two roots of the determinant A. - If therefore in this 
passage from — 00 to zero, as many variations of sign are lost 
as is indicated by the highest power of m^, the values of m* fowrvd 
from the determinantal equation must he all real and negative. 
It wiU also follow that the roots of each of the series of determi- 
nants are alj real and negative, and that the roots of each separate 
the roots of the determinant next above it. If, however, the 
proper number of variations of signs be not lost in the passage 
from m* = — 00 to m* = 0, it does not follow that the values of m' 
are not real and negative. 

4. If the proper number of variations of sign has not been 
lost in this passage from m* = — 00 to m' = 0, this proposition 
does not leave us without information as to the nature of the roots. 
We infer that the number of real negative values of m* is equal 
to or exceeds the number of variations of sign lost by an even 
number and unless the number of variations be even the system 
is unstable. 
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If we do not object to the labour of expanding the determi- 
nants, we might extend this theorem to determine the positions 
of the negative values of tr? as well as their number. The number 
of real negative values of m* between m' = — a and m'* = — ^ is 
equal to, or exceeds by an even number, the number of variations 
of sign lost in the series of determinants. In this form the 
theorem resembles Fourier's theorem in the Theory of Equations. 

5. The converse of this proposition has not been proved. 
If in the passage from m* = — oo to m* = 0, the number of varia-» 
tions of signs is unaltered, it is not true that the values of m* 
cannot be real and negative. Thus in the simple case 






= 0, 



where a,,, a,^^a^ are aU positive, no variations are lost, yet if 
^12 > ^^11 + ^^« *^® values of m* are real and negative. And if 
^18 = ^'^ + ^^^ ^^ roots are equal and negative. It will be 
noticed that the minors in this last case are not zero. 

6. In order to discover the meaning of these losses or gains 
of changes of sign, it will be convenient to make such changes of 
the co-ordinates as will simplify the dynamical determinant as 
much as possible. Let us write 

If we now change our co-ordinates by writing for 6, (f>, &c. 
linear expressions of some new co-ordinates, we know that we can 
clear this expression of all the terms containing the products. 
We know also that this can be done in an infinite number of 
ways. We thus have 



XT ^1 ^1 

»^0 ■" ^11 2" "^ ^22 "2" "^ • • • 

where the symbols 0^^ (f>^, &c. represent the new co-ordinates. 
Again, let us consider the expression 

the coefficients are here the values of P, Q, &c. in the general 
expression for T, Art. 1, Chap. IV. when 0, ^, &c. are all put zero. 
But since T is necessarily positive for all values of f^, 17'. ..f, ^,... 
it follows that W^ is positive for all values of ff, ^\.. Hence 
by a well-known theorem, we may by a real linear transform- 
ation of the variables clear the expression W^ also of the terms 
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containing the products fl', <f>, &c. and can make the coefficients 
of the squares any positive constants we may please. We thus 
have 

^^^ 2 2 



It may be shown that we cannot in general clear the expres- 
sion for T of all the terms containing 0<f)', 0'(f>, &c. unless (7^,— ^^81=0, 
&c. by substituting for 0, 0, &c., any linear functions of other 
variables. As this is only a negative result of which no further 
use will here be made^ it is unnecessary to supply the demon- 
stration. 

7. Using these simplifications the determinant A will now 
take the simpler form, 



A = 






= 0, 



^here o^a = (7^ - C^j, &c. 

K we form from A the series of subsidiary determinants 

A, Aj, A,..., 

terminating with any positive constant, we see that when m' = — oo, 
these subsidiary determinants are alternately positive and negative, 
and when m* = 0, they become 

which are all positive if o^, a^... are all positive. Hence if 
ttu, a^... are all positive, the proper number of changes of sign has 
been lost, and therefore the roots of the dynamical equation A = 
satisfy the condition of stability. 

If G^> fl^M> &c. are all positive, we see that V^ is a minimum for 
all variations of 0^, <f>^, &c. and therefore for all variations of the 
original co-ordinates. If a^^ a^, &c. are not all positive, there will 
be as many variations of sign lost as there are positive quantities 
in the series O-^, a^, &c. In this case V^ is a minimum for some 
variations of 0, <f>, &c. and not for others. It is shown in the 
appendix to Williamson's Differential Calculus^ that n independent 
conditions are necessary that a quadratic expression of n variables 
should be always positive. These are given in the form of deter- 
minants and may be briefly summed up, in the statement that 
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= 0, 



is always positive if the roots of 








-^11 + X, -^12 > 


-^I8» 


&C. 




-^i2> -^sa "^ ^> 


-^28» 


&c. 




-^18 » -^28 > 


-488 + >'» 


&c. 




&C. &C. 


&C. 


&c. 



are all real and negative. 

8. We may now put the proposition of Art. 3 into another 
form. Let L he the general expression for the excess of the kinetic 
energy over the potential energy of a dynamical system in term^ of 
its n co-ordinates x, y, &c. Let this system be moving in steady 

dx dv 
motion with constant values o/ jt > in" > &c. Then if hbe a maxi- 
mum for all variations ofx, y, &c. keeping jt > ^Sx > &c. unchanged, 

then that steady motion is stable. 

If however only t of the n conditions necessary to mxike L a 
wuodmum be satisfied, then the number of roots of the dynamical 
equation which satisfy the conditions of stability is equal to v or 
exceeds r by an even number. There cannot be stability mdess n — r 
is an even number. 

If the system be oscillating about a position of equilibrium, 

-jiy-r: y &c. are all zero, and this leads at once to the condition, 
dt dt 

that the equilibrium is stable if the potential energy is a minimum. 

9. In this reasoning, we have for convenience excepted the 
case in which two successive determinants in the series A, A^, A^... 
vanish for the same value of m*. But this exception is of no real 
importance, for we may change these determinants into others 
whose constituents are very slightly different from those of the 
given determinants but which are such that no successive two of 
the series have a common root. In the limit, therefore, when 
these arbitrary changes of the constituents are indefinitely small, 
the roots of the series of determinants will still be real under the 
same circumstances as before, and the roots of each will separate, 
or coincide tuiih, the roots of the next above it in the series. 

To show that these changes are possible, let us consider the 
row of determinants beginning at the last. The determinant A„ is 
a positive constant, the next A„_j is m* + a^^. Proceeding thus, 
suppose we arrive at two determinants which we may call Aj and A 
which have a common root. If we now change the constituents 
Ou, a,3, ttjg, &c. into a^ + So^, a^^ 4 Sa^^, &c. we do not alter Aj, but, 
except for the root m* = 0, we do alter A in an arbitrary manner. 
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When for example a, is altered, we alter a constituent both in the 
first row and in the first column. Since 

AA' = AA -A A , 

where A' is the determinant formed from A by omitting the first 
and third rows and columns, we see that when A and Aj both 
vanish, the product AjjAgi and therefore by Art. 2 both Ajg and 
A31 must vanish. The determinant A is therefore altered by 
— Aj (Sajg)^ which does not vanish, since A, is by hypothesis finite 
for the particular value of m' under consideration. 

If any determinant of the series vanishes when m* = 0, it is 
clear that one of the quantities a^, a^... must be zero. If we 
replace this by any small positive quantity, the argument will 
apply as befora 

10. It is important to consider the effect of equal roots on the 
test of stability given in Art. 8 of this chapter. 

In this case we know that the roots of the minor A^ separate 
the roots of A. If therefore A have two negative equal roots, it is 
clear that Au must have one of them. In the same way A^, A^, 
&c. must also all vanish for this value of m. Since 

AA' = A A -A A , 

it follows as in Art. 2 if A and A^ both vanish, that Ap^ and A^^ 
also vanish. Hence all the first minors of the determinant A 
vanish. This is the case considered in Art. 5 of Chap. I. The 
equal roots instead of introducing into the expressions for the co- 
ordinates terms which contain ^ as a factor merely' render two of 
the coefficients, instead of one, indeterminate. 

In the same way if the equation A = is satisfied by three 
values of m^ equal to the same negative quantity, the equation 
A^ = must have two of them, and its principal minor must have 
one of them. Reasoning as before we see that all the second 
minors of A must be zero. This is again the test that there should 
be no terms which contain ^ as a factor. 

The presence therefore of equal roots does not in the theorem 
of Art. 8 affect the stability of the motion. 

When a system is disturbed from a position of equilibrium 
whether stable or unstable, the roots of the fundamental determi- 
nant are separated by its minors in the manner described in 
Art. 8 of this Chapter. By what has just been proved, we see 
that if the fundamental determinant have equal roots, whether 
positive or negative, these do not introduce into the integrals 
teims which contain < as a factor. 
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11. The proposition in the last article may be made more 
general. If the fundamental determinant be reduced to the form 
indicated in Art. 7, we shall show that if A vanish for two equal 
negative values of vfC which are numerically greater than the 
greatest negative quantity in the series o^j, a^y a^^ &c., then 
these equal roots will not introduce any terms mto the solution 
with ^ as a factor. If a^^, a^, &c. are all positive, this reduces to 
the propositi<»x prov'ed in the last article. 

Following the same notation as before, we have 

If neither Ajj nor A,, are zero, they must have the same sign 
when A vanishes for a negative value of m\ For their product 
is equal to A^^Ajj which has been proved in Art. 2 to be positive. 
Hence all the leading first minors, viz. Ajj, , A^, &c. must have the 
same sign for any negative value of m^ which makes A vanish. 

By differentiation we have 

^ = 2otA,j + a„A„ + ... - a,^„ f 2mA„ + ... - &c. 

But we have also 

A = (w* + aj Aji + aj^m A,, + ... 
A = - a^^m A^ + {m* + a^ ^m + ••• 
&c. = &c. 

Hence if n be the highest power of m occurring in A, we have 
»» ^ = nA + (♦»* - «m) K + («*" - « J K - &c- 

J A 

If then A and -j— both vanish for any negative value of tn" 

greater than the greatest negative quantity in the series a^^, a,,, 
&c., we have the sum of *a number of quantities all of the same 
sign equal to zero. This requires that each should be zero. We 
have therefore \^^0, Ajg=0, &c. The rest of the proof is the 
same as before. 

J A 

By differentiating the expression for -^ and substituting for 

-^* , -j-^*> &c. their values in terms of their first leading minors, 

we may extend this proposition to the case in which the funda* 
mental determinant has three equal roots, and so on. 
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12. A sphere is suspended hy a string OA from a fixed point 
O, and is set rotating about a vertical diameter which is in the 
same straight line as the string with an angular velocity n. A 
small disturbance is given^ determine if the steady motion is 
stahle. 

Let be the origin, and let the axis of z be vertically down- 
wards, let Ix, ly, I be the co-ordinates of A, the point at which the 
string is attached. Let G be the centre, and let af, arj, a be the 
co-ordinates of G relatiye to A. Then, exactly as in Chap. IV. Art. 
22, the Lagrangian function may be shown to be 






x-^^%'+r^v] 



the mass being taken unity. 

Putting x' = ^> r> Vf ^\ y* ^ zero, we see that Z is a maxi- 
mum when X, y, f, 17 are zero, the steady motion is therefore stable 
for all values rf w. 

If we put i* = I a*, and m* = — \', so that x, y, &c. are all repre- 
sented by terms of the form 2-4 cos (\^ + a), we may, by the 
methods of the last chapter, prove 



(^•-De*-*-!)-!^- 



This equation, whatever may be the sign of w, has two positive 
and two negative roots. All four give stable oscillations. . 



R.A. 
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If tJiA energy of the system he a maodmum or minimum under certain 
conditions, the motion whetJver steudy or not is stable. Arts. 1 — 3. 

When the motion is steady^ it wiU he also stable if a certain function of 
the co-ordinates called Y + W is a minimum. Art, 4. 

If there he only one co-ordinate which enters into the Lagrangianfunctionf 
except as a differential coefficient, this condition is necessary and 
sufficient. Arts. 5, 6. 

Additional conditions when there a/re two co-ordinates. Art, 8. 

1. When a system is oscillating about a position of equilibrium, 
it is well known that we may determine the stability or insta- 
bility of the equilibrium by what we may call the " energy crite- 
rion." This criterion may also be sometimes used when the 
system is oscillating about a state of steady motion. 

Let E be the sum of the kinetic and potential energies of the 
system. Then throughout any motion of the system we have 

E=-h, 

where A is a constant depending on the initial conditions. If 0, <f>, 
&c. are the co-ordinates of the system, E is sl known function of 
0, ff, ^, <f)',&c. Suppose that some of the other first integrals of the 
equations of motion are known. Let these be 

F, (0, ff, </>, <^', &c.) = C, 

K (^. ^> ^> ^^ &c. 

«&c. 

the time t being absent. 

For the purposes of this proposition let us suppose 0, ff, <f>, <f>\ 
&c. to be separate variables unconnected with each other except 
by the equations just written down. 

If E be aii absolute maximum or an absolute minimum for 
all variations of 0, ff, &c., those corresponding to the given motion 
making E constant, then that motion is stable for all displacements 
which do not alter the constants (7^, C\, &c. 
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If this proposition be not evident, it may be proved by elimi- 
nating as many of the letters as possible. If 0, 6\ &c. be the re- 
maining co-ordinates we have 

Let h be the value of E in the given motion, and let the 
system be started in some slightly different manner so that 

E=-h + U. 

If JF be a maximum along the given motion, then any change 
whatever in ^, 6\ &c. decreases E. Hence ^, 0\ &c. cannot deviate 
so much from their values along the given motion that the change 
in E becomes greater than hh. 

2. Let us apply this principle to a system of bodies which 
moves in steady motion with some co-ordinates 0y <^, &c. such that 
their differential coefficients ^, <f>\ &c. are constant, and the re- 
maining co-ordinates f, ^, &e. themselves constant. Let us further 
suppose that the energy is a function of ffy <f>', &c., but not of 
0, <!>. By Lagrange's Equatioiis we have the integrals 

dT^ ^ dT ^ ^ o 

d&^ ^' d^'""^2><»c. 

It is clear that the system can describe any one of a number 
of steady motions, which we have already called parallel motions, 
and which are determined by 

^=;?, <^' = j, &c. 

f = flf, i7 = /3, &c. 

4 

where /?, g, &c* a, ^, &a are constants which satisfy all Lagrange's 
equations. 

We have thus as many relations between these constants as 
there are co-ordinates ^, 77, &c. Let the system be started with 
any initial conditions we please, then the constants (7^, CL, ... are 
given. These being known we have as many relations between 
the constants of steady motion as there are co-ordinates 0y ^, &c. 
The steady motion is therefore determined. If the energy of this 
initial motion is nearly equal to that of this steady motion, and 
if it be a maximum or minimum as explained above, then the 
system will never deviate far from its corresponding position in 
the steady motion, and this steady motion may be called stable. 

3. Example. A top is set spinning on its point on a perfectly 
rough horizontal ground, with its aods inclined to the vertical, find 
the condition of stability. 

Let be the inclination of the axis OG of the top to the 
vertical OZ; yjr the inclination of the plane ZOO to a vertical 

G— 2 
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plane fixed in space, and <f> the inclination to a plane through OG 
fixed in the body. Let be the apex, the centre of gravity 
which lies in OG,h= 00. Let -4, il, (7 be the principal mo- 
ments of inertia at 0, and if the mass of the top. We have then 

^=f (^' + t'«>s^" + ~(^» + sin«^f'«)+%Acos^. 

By Lagrange's equations we have the integrals 

(j) + '^' cos ^ = n, 
CVi cos ^ + -4 sin'0^' = m, 

where n and m are two constants, the former representing the 
angular velocity of the top about its axis, and the latter the 
angular momentum about the vertical. If we now eliminate 
^' and '^y we find that J? is a minimum when = a,i£ 

G^rf > 4iMghA cos a, 

which is the result given by other methods, 

4. The theorem of Art. 2 may be put into another form. 
Let the kinetic energy be 

Then since 6, if>, &c. are absent from the coefficients, we have 
the integrals 

T^^ + r^f + ... = c^ - T^^ - r^v - &c. 

&c. = &c. 

For the sake of brevity let us call the right-hand sides of these 
equations G^ — X, C, — F, &c. Since T is a homogeneous func- 
tion of ff, <f>\ &a, we have, as in Chap. iv. Art. 21, 

+ K(c^i+^)+H'(Ci+r)-h... 

If we substitute in the second line the values of ff^ <]>, &c. 
found by solving the linear equations just written down, we have 
the determinant 

0, (7, + Z, G,+ Y, &c. -^ 
0,-X, Tee, Te^. &c. 
G^— Y, Tii, Ta^^ &c. 



&c. 



&c. 



&c. &c. 
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where A is the discriminant of the terms in T which contain 6\ (f>\ 
&c. This determinant is unaltered by changing the signs of X, Y, 
&c. and is a quadratic function of C^, (X, &c., X, Y, &c. Hence 
the tenns C^X, G^Y, &c. do not occur. If then we put 



F=- 





0.. 






&c. &c. 



0,. 

&c. 



&c. 
&c. 
&c. 
&c. 



2A 



we have 



T=B^X + B^^rf' +&C.+ W, 



where B^^, &o. are independent of C^,, 0^, &c. Now Tis essentially 
positive for all values of the variables, and therefore for such as 
make (7^, C^, &c. all zero. Hence the terms involving f, r}\ &c. 
are together a minimum when ^, rf', &c. are all zero. The coeffi- 
cients B^^, &c. may all be treated as constants since f ', rj, &c. are 
all small quantities. 

If F be the potential energy, we have therefore the following 
rule. J/* W + V be a minimum for all variations of f , rf, &c. tiien 
the steady motion is certainly stable. It should be noticed that 
TF+ F is a function only of f, rj, &c. the co.-ordinates which are 
constant in the steady motioa. 

5. If the energy be a fimction of one only of the co-ordinates, 
though the differential coefficients of all the others enter into its 
value, this condition is sufficient and necessary. 



Let ^ be this co-ordinate. Then by vis viva we have 
Differentiating we have 



BiS^r + W+V=k 



^i^^ + — 5I — ="^- 

This equation must be satisfied by the steady motion repre- 

sented by ? = a. The second term ^ ,^ must therefore 

vanish when f = a, so that W+Vis a maximum or minimum. 
To find the oscillation let us put ^= a + a?, we find 



-B«|- 



H- 



(TF+F) 



d^ 



1ar = 0, 
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where the square bracket implies that f is to be put equal to a 
after differentiation. By the same reasoning as before B^^ is 
necessarily positive, and the motion will he stable or umtahle oc- 
cording as (W + Y) is a minimum or mxiocimum. 

6. If we refer to Art. 21 of Chap. iv. we see that this function 
Tr+ Fis the value of the modified Lagrangian function L' when 
f , 17', &c. are all put zero, and the sign of the whole function 
changed. It therefore follows by Chap. V. that when TF4- F is a 
minimum the steady motion is stable. The "energy criterion of 
stability," as far as it applies to steady motion, may therefore be 
deduced from that given in Chap. v. Art. 8. But the mode of 
demonstration adopted in that chapter gives us more information 
as the nature of the motion, while the modes of application to 
examples of the two criteria are quite different. The energy 
criterion may also be sometimes applied to determine the stability 
of a motion which is not steady. 

[The relation between the theorem in this Chapter in which 
E= T-\- V is made a minimum to that given in Chapters rv. and V. 
in which L = T— V is made a maximum may be more distinctly 
perceived by the following statement. 

Let Xy y, &c. be the co-ordinates of the system, and let L be 
the Lagrangian function* so that L = T— F, then by Art. 3 of 
Chap. IV. the co-ordinates in steady motion satisfy the equations 

Here L is expressed as a function of x, y, &c. x, y\ &a 

Suppose some of the co-ordinates as 0, <f>, &c. are absent from 
the expression for L, so that i is a function of ff, <f>\ &c., the re- 
maining co-ordinates, viz. f , 77, &c. and their differential coefficients. 
Then if we form the modified Lagrangian function as in Art. 21 
of Chap. IV. the equations (1) of steady motion become 

f-". f-»'^- » 

Here, as in the Hamiltonian equations Cj, C^, &c. are the 
d, (f>, &c. components of momentum, and L' is expressed as a 
function of C^, C^, &c. f, 97, &c. |', rj\ &c. 

dT 

* Let M, V, &c. be the ar, y, &c. components of momentum, so that m=3-,, &c. 

ax 

and let ff be the Hamiltonian function. Then H^T-^ V and we easily deduce 

from the Hamiltonian equations that, in steady motion, 

dll ^ dH ^ ^ 
-7- = ", - ,- =0, &c. 
dx dy 

Here H is expressed as a function oix, y, &c. w, r, <tc. 
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But i' = r-C/-(7,f-&c.-F 

by Euler's theorem on homogeneous functiona In steady motion 
f , 17', &c. all vanish, hence the equations (2) become 

^dJW+Jl d{W+V) _ 

d^ "' dfl -"'®^ w. 

where W is the value of T when f , rf\ &c. are all put equal to 
zero. Here TFis expressed as a function of C^, 0,, &c. f, iy, &c. 

It is shown in Chap. v. that if the Lagrangian function ex- 
pressed as required in equations (1) be a maximum the motion 
is stable. It is shown in this Chapter that if the function W+ V 
be a minimum the motion is stable.] 

7. To find ihe condition of stability when the Lagrangian 
function is a function of two only of the co-ordinates, though the 
differential coefficients of all the others enter into its value. 

Let f, iy be these two co-ordinates, then the modified Lagran- 
gian function as explained in Art. 20 of Chap. iv. will be a fiinc- 
tion of f, iy, t , 17 only. 

Let the steady motion be given by f = a, 17= A with the 
corresponding values of the other co-ordinates 0, <f>, &c. Then a 
and /8 are constants. Let f = a + a?, i7 = )8 -f-y, and let us expand 
the modified Lagrangian function in powers ^ x, y. Neglecting 
the terms of the first order, as they only give the steady motion, 
let 

-^ G^xx -h G^y' + G^cd + G^y\ 

Also let E=^ Cfl— Cia- Then the condition of stability is that 
the roots of the following equation should be of the form /8 ^/— 1. 



B^w^ — A^ — Em, B^ni? — A^ 



= 0, 



If A'^B^B^'-BJ 

A 



^B^B^-B^ \ 
and = Jui?„ -f- A^y^ - 2^i,Bia 
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be the two discriminants and the oth^r invariant, this leads at 
once to the conclusion that the motion is stable only when 

(1) A is positive, 

(2) JST* - is positive and > 2 VAA'. 

If Ayj^ + A^y + -4„ |- is a maximum when x and y are zero, 
the two conditions are obviously satisfied. 

This condition may be otherwise expressed ; if L' be the modi- 
fied Lagrangian function, a steady motion is given by 

%'"■ f -». ?=»■ '•-• 

This motion will be stable if for the values of f, 17 thus found, 

^^^ [d^dv' d^dv) |df dv" "^ dri' d^"" d^drid^'dr{] 
is positive and greater than 



2 



{d^ud'^L d^LM {d^rd^n dT) 



yf drf d^dri] '\dl^^ dri\ d^'d^} * 



i 



8. The nature of the motion when thus reduced to depend on 
two co-ordinates may be illustrated by geometrical reasoning. Let 
the position be defined by two co-ordinates x, y which are zero 
along the steady motion. Let these be regarded as the co- 
ordinates of a point P referred to any axes. Ihen the motion of 
P exactly represents that of the system. 

Let us construct the conies 

2 2 

A ^ A , A V 

^U2+A2^ + ^a2|- = a, 

8 2 

-Bu|+5i2''^y + £22| = J. 

Then, exactly as in Arts. 11 — 14 of Chap. IV., it will be found 
convenient to transform the co-ordinates by writing 

If ft be the modulus of transformation we have ft = a^^ — a}>^. 
It is easy to see by actual substitution, if E' be the difierence of 
the coefficients of j)^' and pq^ that 

E = iiE. 
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If the transformation be from one set of oblique co-ordinates 
to another, let o), <o be the angles between the axes. We then 
have 

E' ^ E 

sin w' sin 0) * 

Transforming the axes to the common conjugate diameters^ the 
conies become 

A '^4-A *t^a 

the signs of a and h being such as to make these conies real. The 
equation to find m becomes 

(£>« - A^) {B^m' - ^ J + E^'m' = 0. 

It is therefore necessary for stability that the conic a should be an 
ellipse as well as the conic 5, It is also necessary that. 






both roots having the same sign and the inequality being nu- 
merical. 

Let OP, OP ; OQy OQ he the common conjugates of the two 
conies, this condition then becomes 

E area of conic 6 OQ OQ 



sin©* ^v^ ^ OP^ OF' 

When the system describes an oscillation with one period, i,e, 
an harmonic oscillation, the path of the representative particle is 
easily seen to be 

( BJW - A,;')f+ (5/m« - A^') ^ = constant. 

The harmonic paths are therefore ellipses. It also appears that 
the two ellipses which represent the two harmonic vibrations and 
the two ellipses a and 6 have, all four, a common set of conjugate 
diameters. * 



CHAPTER VIL 



Any small term of a high order, {fits period is nearly the same as ihat of 
an oscillation of the system, may prod^ice important effects on the 
m>agnitude of the oscillation. Art, 1. 

Origin of such terms, with an example. Arts, 2 — 3. 

Supposing ths roots of the determinantal equation to satisfy the conditions 
of stability to a first approximation, yet if a commensurable relation 
hold between these roots it is necessary to examine certain terms of the 
higher orders to determine whether tJiey will lUtimately destroy the 
stability of the system. Art, 4. 

If a certain relation hold among the coefficients of these terms, they will 
not affect the stability of the system, but only slightly alter the periods 
of oscillation. Arts, 5 — 7. 

Examples, the fi/rst taken, from La>grange^s method of fimding the osciUa* 
tions about a position of equilibrium. Arts, 8—9. 

If the coefficients of the equations of motion slwuld not be strictly con- 
stant, hU only nearly so, the stability will not be affected, unless the 
reciprocals of their periods liave commensurable relations toith the 
reciprocals of the periods of oscillation of the system. Art, 10. 

1. If we understand that a motion is called stable when any 
small disturbance does not cause the system to deviate far from its 
undisturbed motion, it is clear that we cannot be certain of the 
stability without examining the terms of the second order. It is 
possible that some of these may have their periods so timed that 
their effects may accumulate until the motion is changed. 

Returning to the equations of Art. 3, Chap. I. we shall have on 
the right hand, instead of zero, a series of small terms of orders 
higher than the first. To find a second approximation, we substi- 
tute the values of a:, y, &c. given at the end of Art. 3, in these 
terms. 

They will therefore take the form Ne***'Sind will produce in 

or, y, &c. terms of the form -^y-^ e***, where W is of the same order 

at least as the term considered, and/(m) Jias the same meaning as 
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in Chap. I. These will have to be expressed in trigonometrical 
real forms, but it is unnecessary to exhibit the process, for we see 
at once that no small term or force (whatever it may be called) of 
a high order can aflfect the stability of the motion unless it makes 
f{n) very nearly or exactly equal to zero. In this case its period 
is very nearly or exactly equal to one of the periods of the motion 
given by taking terms of the first order only. 

A remarkable use of this principle was made by Captain Kater 
in his experiments on the magnitude of gravity. It was important 
to determine if the support of his pendulum was perfectly firm. 
He tells us that he had recourse to a delicate and simple instru- 
ment the sensibility of which was so great that had the slightest 
motion taken place in the support it must have been instantly 
detected. The instrument consists of a steel wire the lower part 
of which-, inserted in the piece of brass which serves as its support, 
is flattened so as to form a delicate spring. On the wire a small 
weight slides by means of which it may be made to vibrate in the 
same time as the pendulum to which it is to be applied as a test. 
When thus adjusted it is placed in the material to which the 
pendulum is attached, and should this not be perfectly firm its 
motion will be communicated to the wire, which in a little time 
will accompany the pendulum in its vibrations. This ingenious 
contrivance appeared fully adequate to the purpose for which it 
was employed, and afforded a satisfactory proof of the stability of 
the point of suspension. See Phil. Trans., 1818. 

2. Since the term Ne"** is obtained by compounding the 
different terms in the values of a?, y, &c. it is clear that 



• • • 



n =pmj^ + qm^ 

where p, q, &c. are positive integers whose sum is the order of the 
term. It is therefore only when the roots of the dynamical equa- 
tion/ (in) = are such that a linear relation of the form 

pm^ + qm^ + ... = m^ very nearly 

exists between them, that we may expect to find important 
terms among the higher orders. The order of the terms to be 
examined will he p + q + ..., and unless this be also small, the terms 
will probably remain insignificant. If the root m^ should occur 
twice in / (m) = it is clear that the divisor / (n) will be a small 
quantity of the second order, and the term may be said (as in 
the Lunar Theory) to rise two orders. 

3. To take an example, let us suppose a particle to be describ- 
ing an ellipse about a fixed centre of force in one focus. If dis- 
turbed it will describe a slightly different ellipse. Since 

r = a{l — ecos (n^ + € — ft)) + ...}, 
d=.nt+ € + 2e sin {nt-\- € — (o) -\' ... 
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we see that a slight change in the elements will cause variations 

in r and of the period — , an additional variation in of the 

form iSn + he and an additional variation in r of the form Sa« All 
these variations should by Art. 3 of Chap. i. be indicated by 
expressions of the form 

where the values of m are the roots of the equation f{fn) = 0. 
The roots therefore of the equation f{m) = for Sr are m = 

and ± nV — 1, and for 85 are w = 0, and ± » V — 1. We there- 
fore infer that any small disturbing causes of the second order 
whose periods are nearly equal to that of the particle, will cause 
important inequalities in both Sr and hO, and (since f{m) =Ohas 
two roots equal to zero) any term of long period will rise two 
orders in h0. 

4. If the roots of the subsidiary equation are such that the 
relation 

holds accurately, the solution changes its character. We have 
now in the value of a; a term of the form te'^^^. Unless the real 
part of m^ is negative, this indicates that the system will depart 
widely from the motion which we took as a first approximation. 
We must therefore modify our first approximation (as in the Lunar 
Theoiy) by including in it the terms which produced these im- 
portant effects. We may then enquire how far this modified first 
approximation indicates that the motion is stable or imstable. 
When these terms are included the equations to be solved are in 
general no longer linear, and it is sometimes impossible to find a 
solution sufficiently accurate to serve as a first approximation 
throughout the whole motion, 

5. In some cases, however, the oscillations may still be 
represented by expressions of the form 

&c., 

where the values of n^^ n^.. differ but slightly from the roots of 
the equation f{rn) = 0. Let us investigate the condition that this 
should be true, and also determine whether the changes in the 
values of m^, m^. . . are sufficient to affect the stability. 

Suppose that we have completed our first approximation, and 
find on proceeding to a higher approximation that the terms 
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present themselves on the right-hand side of the first set of 
equations in Art. 3, Chap. I. These terms are supposed to have 
arisen from several relations of the form 

If these terms can be included in the first approximation by 
writing n^y n^, &c. for m^, w^, &c. we have, by substitution in the 
differential equations, certain equations connecting n, M, M\ &c., 
whose left-hand sides are the same as those used in Art. 3 with 
Wi written for wij, but on the right-hand sides we have instead of 
zero the quantities ^, jY/, &c. The test of the success of the 
process is that these modifications in the values of m must satisfy 
the same relations as before. 

Now N^ iV/, &c. are all at least of the second order of small 
quantities, hence up to terms of the first order the ratios M^y Jlf/, 
&c. will be the same as before, so that we may put 

following the same notation as in Art. 3, Chap, i. We also have 
Let Wj = TWj + Swi we find 

Similarly 

&c. = &c. 

It is evident, by the theorem of determinants alluded to in 
Art. 3, that these are symmetrical expressions. 

6. We may conveniently express these results in the form of a 
rule. 

Suppose we have to a first approximation 

Eliminate from the diflferential equations all the variables except 
a? in the usual manner. This may be done by performing on the 
several equations the operations represented by the minors a, a, 

&c., -J- being written for m. Let the equation thus found be 
(It 



fQ-^F^'-' + P^^''- 
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Then all these terms can be included in the first approximation, 
provided 

satisfy the relations 

pSm^ + qBm^ + . .. = Sm^, 

&c. = &c. 
which exist among the roots of the dynamical equation. 

7. The general results we have arrived at may be summed 
up as follows. Though some of the terms of the higher orders 
may affect the. magnitude of the oscillation, yet no term will arise 
to aflfect the stability of the motion unless there be some relations 
between the roots of the dynamical equation of the form 

where p, q, &c. are all integers. Even if such relations occur, the 
lowest order of the term is p + q-^- ...,-and if this be considerable 
the term will not produce any important effects until a con- 
siderable time has elapsed. If a certain relation, just found, hold 
among these terms, their only effect is slightly to modify the 
periods of oscillation, without altering the type of niotion. 

8. As an example, let us consider a system of bodies to be 
Oscillating about a position of equilibrium. We know by Lagrange's 
general solution, that the equation f{m) = is of an even order. 
Its roots are of the form 

77ij = aV— 1, mj=: — aV— 1, m^=^^^—l, &c. 

Whatever the numerical values of these may be we have 

m^^-m^-Vm^^m^y m^+m^ + m^^m^, &c. 

so that the small terms of the third, fifth, &c. orders might affect 
the stability of the oscillation. But we shall now show that they 
only affect the periods of oscillation, and not the stability of the 
system. 

Since both sides of Lagrange's equations must be of — 2 di- 
mensions in time and the impressed forces are also of — 2 dimen- 
sions, it is clear that these terms must consist of powers of x, y, &c. 

~j^9 "j^f <^c* ^^^ products of an even number of factors of 

■-tt , -^., &c. We know also, by Lagrange's solution, that the 

co-ordinates take the form 

X— M^ cos (cLt + XJ + &c. 
y = Ml cos {it + \) + &c. 
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The minors a, a', &c. are also all even functions of -^ , hence 
the equation found after elimination is of the form 

P cos {oi + \) + ^cos 03« + aa) + ... 



/(!)'= 



Replacing aV— 1, *-aV^, &c. by Wj, m^, &c. we find by 
"Art. 6, 

P P 

Since /'(m) is of odd dimensions, and 7Wj=— w, we clearly have 
Bm^=—&m^ and therefore the test is satisfied. 

9. [As another example let us apply the rule of Art. 6 to 
some very simple case which will involve no algebraical substi* 
tutions of any length. 

The motion of a simple pendulum under the action of gravity 
may be made to depend on the equation 

g + a'x=^x» (1), 

where a and 13 are two constants and x is the inclination of the 
pendulum to the vertical which is supposed to be small The first 
approximation to the motion is 

a = if jfi^i* + Jif^e"*** (2), 

where m^ and m^ are the roots of the equation m* + a' = 0. Our 
object is to ascertain by help of the rule given in Art. 6 whether 
the small force represented by Bx^ renders this first approximation 
unstable or merely slightly alters the numerical values of m^ and 



Wlj. 



The two roots are connected by the relation 

*mj + «2a = (3). 

Substituting the value of x on the right-hand side, we have 
Hence by the rule in Art. 6 

These clearly satisfy the relatioa 

Sm, + Sotj = 0, 
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and therefore the first approximation taken above, so far as the 
disturbing force Bo^ is concerned, is stable. 

If the small force had been ^\-rA instead of ^sf^ it is easy 
to see in the same way that 

_ 3/3Jlf/Jtf.m>, _ Zm,M:m,fn: 

^> — ^,M, ' *"« — 2^;]^ — ' 

so that the relation Srn^ + ^m^ = would not have been satisfied. 
The first approximation taken above is therefore not sufl&ciently 
accurate to serve as a first approximation throughout the motion. 

In this example we have considered the effect of a small force 
of the third order in disturbing the stability of the motion given 
by equation (1). The same equation will obviously occur in many 
other cases of motion. For example, let a particle describe a cir- 
cular orbit about a centre of force situated in the centre. If 
slightly disturbed the equation giving the disturbance x in the 
radius vector takes the form 

where a, /3 and 7 are constants. Similar remarks will therefore 
apply to this case also.] 

10. When the coefficients of the equation of motion are not 
strictly constant, but yet do not vary much, then we may transpose 
the small variable parts of these terms to the right-hand side of 
the equation, and treat their products by the differential coeffi- 
cients of the co-ordinates as small quantities of the second order. 
Suppose the variable part of one of these coefficients to be p sin nt, 
where p is small, and let/(m)=0 be the equation givmg the 
periods of oscillation of the system when the coefficients are taken 
constant. Then it is clear that unless n is nearly equal to the 
sum or difference of two values of m, this term cannot rise into 
importance. On proceeding to higher orders we see that these 
terms cannot produce important effects unless some commen- 
surable relation between m and the roots of the equation /(m) = 
should be very nearly satisfied. 

11. It should be remarked that when the coefficients are not 
constant it is not a sufficient test of stability that they should 
always satisfy the conditions of stability obtained by giving them 
their instantaneous constant values. Thus if the- equation of 
motion were 
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the coeflScient of x is always positive, yet as the equation is satis- 
fied hj x = af^t,x may become as great as we please. 

Even if the coefficients are nearly constant, we must yet ex- 
amine, by the rules just given, if their small changes are so timed 
as gradually to increase the oscillation until the divergence from 
the given motion is no longer small. 

[Suppose a system to be oscillating so that its motion is de- 
termined by the equation 

where j is a known function of t, which during the time under 
consideration always lies between ^8* and /8'* the latter being the 
greater. Let the system be started with an initial co-ordinate x^ 
and an initial velocity a?/ in a direction tending to increase x. 
It may be shown that the system will begin to return, i.e. x will 

x^ 
x^ + -^. If 

± m, +m' be two successive maximum values of x, we may also 

show that m cannot be so great as ^.fn, and that the time from 

ODO maximum to the next lies between -5 and ^ ], 



R. A. 



CHAPTER VIII. 



The Hamiltonian Gharacteriatic or PrindpdL functiona wh&n found deter- 
mine at once the motion of the system, from one ffiven position to 
arwiher^ and whether the motion is stoible or unstable. Arts, 1 — 3. 

Examples with a mode of effecting the integration S= iLdt m sm4dl oscUr 
latums. Arts. 4 — 7, 

T/ie Characteristic ftmction supplies the condition t/iat the motion is stable 
as to spojce ordy, while the Principal function gives the conditions that 
it is stable both as to space and time. Art. 8. 

In wliat sense the motion is unstable if eit/ter of the two HamUtoniom 
functions is a minimum. Arts, 9 — 14. 

1. If we had any convenient methods of finding the Hamil- 
tonian Characteristic or Principal function, we might detennine 
without difl&culty the conditions of stability of a dynamical system 
at the same time that we deduce the integrals of the equations of 
motion. But it is very diflScult to discover either of these func- 
tions by an d priori method. We have indeed differential equa- 
tions which they must satisfy, and Jacobi has taught us what kind 
of solution will serve our purpose. But the difficulty of finding 
these solutions is as great as that of solving the equations of 
motion. For these reasons it does not seem necessary to dwell on 
the uses of these functions, 

2. Suppose the Principal function /S' of a dynamical system to 
have been found in terms of the initial co-ordinates 0^ <f> and the 
co-ordinates 6, ^ and the time t. Let the semi- vis viva be given 

by 

where P, Q, jR are known functions of 6, (p. Let P^, Q^, R^ be 
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the values of these when 0^^ ^^ are written for 0, <f>. The final 
integrals of the equations of motion are then given by 

Let the system receive any disturbance at the time ^ = 0, so 
that while starting from the same initial position, its initial ve- 
locities are slightly altered. Let x^'y yl be these initial changes 
of 01 and ^l and let 0-\-x, <f>-\-y be the co-ordinates of the 
system at the time t Then we have 



d0,d0'^ d0,d(f> 
d'S d'S n ^ T> ^ 

^d^0''''dmy^^^^'^^^y^ 

Here x^, yl are small arbitrary quantities, hence x and y will 

be small if none of the ratios of ^^ , -^, -^, or -^ 

dd^d0' d0^d<t>' d<l>^d0' d<f>^d4> 
to the determinant 

^ ^^^^ d^8 
dJ0^d0 d^^d^ dJd^d^ d<^40' ' 
be large. 

If the initial position as well as the initial motion be altered 
we may find, by a precisely similar process, the conditions that 
X and y should be small. If the system have more than two inde- 
pendent motions, we have more than two co-ordinates, but the 
conditions of stability are found in the same way. 

3. If a: and y be small throughout the whole motion from the 
one given position to the other, not only does the system not 
deviate far from its undisturbed course, but the system at any 
instant is also very nearly coincident with its undisturbed place 
at the same time. It is important to notice this, for the word 
*' stability " is sometimes used in a different sense. 

This condition of stability may be put under a form in which 

no reference is made to 8. Let u and v be the components of 

momentum of the system corresponding to the co-ordinates 0,6 

dT dT 

respectively, i,e, let w = ;^/i ^~^» ^^^ ^^^ these be expressed 
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OS functions of 0^, <f>^t 0, (}> and t Then the preceding equatiou 
may be written 

du . dv 

du , dv r> 

where a and fi are small arbitrary quantities. The condition of 
stability is that the values of x and y thus found should be small. 

4. As an example let us consider the case of a projectile. 
If d be the horizontal, and (j) the vertical co-ordinate of the par- 
ticle, we have 

The equations to find z and y are evidently 

1 

— X ^= X 



\y-y:\ 



Hence the systeni continually deviates more and more from its 
undisturbed place. 

6. In order to calculate the form of 8 when a system is oscil- 
lating about a state of motion, it is convenient to choose as co- 
ordinates some small quantities x, y which vanish in the given 
state of motion. Ijet the Lagrangian function be written in the 
form 

where L^ is a homogeneous function of a?, y, a?', y\ Then by a 
theorem of Euler's 

where the 2's imply summation for all co-ordinates. 

As in Art. 9 pf Chap. iv. we have 

d dL^ ^ dL^ 
dt dx ^ dx ' 

a-nd the oscillations are given by 

d dL^ _ dL^ 
dt dx dx " 
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Hence we find 
Integrating we have 

»./..*.. p^]>is [«§.];. 

Thus the integration has been effected, but in order to express 
S as a function of x^, Vo, x, y and ty it will be necessary to find 
x^y y^y X, y in terms of these quantities. 

6. As an example, let the position of the system depend on 
one co-ordinate x and let 

L^L,^-A,x-{'B,x'-^\A,,x^-\'\B^ai^-^C,,xx\ 

where the coefficients are all constants. We then find by the 
process just indicated that -4^ = and 

8^L,t-\-B,{x-x,) + ^na?-^x:) 

mB^^ {a? + x^) { e^ + g"^) - 4;^^ o 

A . • . . 

where vr? = -^ ♦ Applying the criterion of stability we find that 

II 

•j — -r- will finally become small if m is reaL The motion is there- 
dx^dx 

fore unstable or stable according as -d-n, B^^ have the same or 

opposite signs. 

7. If the position of the system depend on two co-ordinates 
X, y, let 

L^L^-^-A.x^-Aji-^B.x'^-B^y' 

+ }^A^o?+A,,xy + ^AJ^-\-^B^x''^B,,x'y' + i,B^y'' 
+ G^xx + C,^y' + C^^yx' + C^y'. 

We then find 
S=Ljt + B,{x''X,) + B,{y''y,) 

^0^,'^•,^^^xy^x,J:)^CJ^^^ 
where 
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If we now express x, y\ x^y yl in terms of a?, y, x^ y^ and f, we 
find for a a fraction whose numerator is a homogeneous quadratic 
function of a?, y, x^ y^,, the coeflScients being Imear functions of 
exponentials of ^, and whose denominator is another linear function 
of the same exponentials. These exponentials become sines and 
cosines when the motion is stable. Thus when the given motion 
is steady the simplest inspection of the form of a will determine 
whether the motion is stable or not. 

Referring the motion to principal co-ordinates for the sake of 
brevity, and writing 2 G = (7^^ — 0,^, we find that a must satisfy 
the differential equation 

This equation is obviously satisfied by such a function as that 
just described. The solution of this equation may be reduced to 
linear equations and thus a may be found. But it is unnecessary 
to dwell on this, for this would be equivalent to returning to the 
Hamiltonian equations* 

8. If we wish to determine the condition that the general 
course of a dynamical system is stable without requiring it should 
be near its undisturbed place at any the same time, it is more 
convenient to use the Characteristic function. Suppose that the 
Characteristic function has in Jacobi's manner been expressed as 
a function of the co-ordinates ^, ^, the constant h of vis viva and 
two arbitrary constants a^, a^. Then 

The relation between 6 and <^, which may be called the equa- 
tion to the path of the system, is given by 

dV , 

where \ is another constant. Let the system be disturbed from 
the same initial position so that the whole energy is unaltered. 
The change in ^ corresponding to any given value of ^ is found 
from 

(PF 
Let A be the initial value of -7—5 , then 

cfF . 
dajd<l> 
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The condition that the path should be stable is that the coeffi- 
cient of Sttj should not be large. 

We might also use the function called Q by Sir W. R. Ha- 
milton, but it seems unnecessary to dwell more on this subject. 

9. The instability of a system may be deduced from the 
Hamiltonian Principal or Characteristic functions, expressed as a 
minimum. Suppose a dynamical system to move from one posi- 
tion A to another J9 in a time ^, then the motion may be found by 

making the first variation of iS= 1 Ldt equal to zero, the time of 

Jo 
transit being constant. The constants of integration are deter- 
mined by the conditions that the co-ordinates have given values 
when t = and t = t To determine whether 8 ia a. maximum or 
minimum or neither we must examine the second variation and 
here we have the assistance of Jacobi's rule. The determination 
of the constants will depend on the solution of equations and may 
lead to several different kinds of motion from A to B. One of 
these will be the actual motion. Let us move B along this until 
one of the other motions coincides or as we may say approaches 
indefinitely near to this actual motion. We have then reached 
a boundary beyond which the integration must not extend if iS is 
to be a maximum or minimum. See Todhunter's History of the 

cPL 
Ccdculus of Variations, page 251. Further -^ , if ^ be a co- 
ordinate, is positive throughout the limits of integration, so that 
8 will be a minimum and not a maximum. 

10. When there are several co-ordinates 0, (f>, &c. which are 
to be found as functions of the time, we may easily show that 
Jacobi's condition is a necessary one, and this is all that we require 
for the next proposition. If the system can move in two ways 
from A to B, then 88=0 along each, and therefore when these 
two are adjacent we have both S/S = and S(8-\-SS) = 0. This 
shows that the second variation can be made to vanish by taking 
one variation through the other. This second variation will then 
be the same as the quadratic term of the series obtained by 
changing the co-ordinates 0, (f> into 4- B0, ^ 4- 8<f>, because we can 
take S^0 = and S*<^ = 0. Hence as the sum af the terms of the 
third order does not, in general, vanish for this displacement, it is 
clear that 8 cannot be either a maximum or a minimum. 

Let the actual motion be from A to By and let a neighbouring 
motion starting from A lead the system to a position G reached in 
the same time along the actual motion before reaching B, Then 
we can show that a variation of the actual motion from A to B 
can be found which makes S^S of any sign. Let P be any position 
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on the neighbouring motion before reaching (7, and Q, one on the 
actual motion after passing C. Then considering Pand Q as fixed 
and also the time of transit, the motion along PCQ cannot 
make S/S= 0; for this condition is known to lead to the ordinary 
dynietmical equations, and it is clear that (impulsive forces being 
set aside) no actual motion can be discontinuous. But there is 
discontinuity at C, for otherwise when the system is started from 
B towards A, two courses would be open to the system on arriving 
at G. Hence the first variations of 8iot an imaginary motion along 
FCQ are not zero, and therefore may have any sign. But since 
the discontinuity at G is of the first order of small quantities, this 
first variation is of the second order. Now the value of S for the 
actual motion is equal to that along the neighbouring motion to 
G and then along the actual motion to B. Hence, P and Q being 
still fixed, variations of the actual motion from A to B can be 
found which make S^8 of any sign, 

11. Let us apply this theory to determine the stability of a 
given state of motion. First let us supposq the given motion to 
be steady and to depend on only two variables. If we use the 
function 8 there will be one co-ordinate and the time, if V two 
co-ordinates. Let the system be disturbed at any moment by an 
alteration of the velocities of its several parts, so that the initial 
position of the disturbed motion is an undisturbed position. If 
the motion be stable the system will oscillate about the un- 
disturbed motion, the oscillation repeating itself at a constant 
interval. It follows therefore by Jacobi's rule that 8 or V cannot 
be a minimum for a period longer than the time of a half-oscilla- 
tion. If therefore /S or F be a minimum for all variations, starting 
from A and ending at JB, where J5 is a position on the steady 
motion reached by the system at an interval as long as we please, 
then the motion is unstable. 

If we give a meaning to the word " stable" somewhat different* 
from its usual signification, we may extend this proposition to 
determine a test of the stability of any motion, whether steady or 
not. All we have assumed is, that, if the motion be not altogether 
unstable, there are some disturbances which will cause the system 
periodically to assume the same positions as it would have done if 
it had been undisturbed, but the interval of these periods may be 
any whatever provided the first be finite. If we use the Character* 
istic function, these disturbances must be such as not to alter the 
constant of vis viva, and if the Principal function, they must be 
such as to bring the system to an undisturbed position in the same 
time. 

* This meaning docs not always agree with the results of Art. 14, Chap^ iv, 
[See also Aits. 17 and 18, Chap, iv.J 
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12. Next let us suppose that as the system proceeds from the 
initial position A along the actual motion, a ceases to be a 
minimum at some position B. The conditions for a minimum 
are df two kinds. Suppose the system to depend on two co- 
ordinates Oy <f>, and let L be the Lagrangian function, then (1) we 

must have -igrj and ^z^a jt?5 "" "^ZTrn ^^^ positive, and (2) it 

must be possible to choose three arbitrary constants which enter 
into a very complicated expression, so that this expression may 
never become infinite between the limits of integration. The first 
condition is clearly always satisfied since the vis viva of any 
system is necessarily positive for all values of d' and <f>\ The 
second condition will fail if there are two neighbouring motions by 
which the system can proceed from A to any position between 
A and B. If this be the mode of failure, it is clear from the 
reasoning of Art. 2 that the conditions of stability are satisfied for 
one kind of disturbance, and that therefore some at least of the 
harmonic motions are stable or oscillatory, though the motion may 
be unstable for a different kind of disturbance. 



13. [These conditions become much simpler when the position 
of the system is determined by one co-ordinate, or when the 
Lagrangian function can be reduced to depend on one co-ordinate. 
Let this co-ordinate be so chosen that it vanishes along the given 
motion, and let us also suppose that both it and its differential 
coefficient with regard to t, are small for all neighbouring con- 
strained motions. Let this co-ordinate be called and let the 
Lagrangian function be 

L=^L, + A,0+ B,ff + 1 A,,0' + 1 BJ' + CJff. 

Then since the Lagrangian equation of motion is satisfied by 
hypothesis when ^ = 0, we have A^ = -B/ where the accent, as 
usual, denotes differentiation with regard to t 

If the system be now conducted from the initial position A to 
any other position B, both on the given motion, by any neighbour- 
ing mode of motion, we have 

5= JL,dt +J(| BJf* + CJff +\aJ^ a 

If ^ = M be any solution of the Lagrangian equation 

J^{BJ^C,JS)^CJ + AJ, 
R. A. 8 
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we may write the function 8 in the form* 

The second term is essentially positive, since B^^ must be 
positive. Hence ^ is a minimum along the given motion unless 

we can so choose the arbitrary displacement d as to make -^ - = 0. 

This gives 0=cu where c is some constant. But must vanish 
at the two limits A and B, hence this choice of is excluded 
unless there is some neighbouring mode of motion by which the 
system could move freely from the given initial position A to the 
position B. The result is that 8 cannot cease to be a minimum 
before the first instant at which some neighbouring motion will 
bring the system (starting from A) into coincidence with some 
contemporaneous position on the given motion. If the given 
motion be steady, it follows that 8 cannot cease to be a minimum 
before a time which is half that of a complete oscillation. 

We thus have a test of stability. If the system depend on one 
co-ordinate and if S be a minimwm when the limits of integration 
are from the initial position A to all positions on the actual motion^ 
that motion is v/nstahle. But if S cease to he a minimum at some 
point C, tlien the actual motion is stable from A to C] 

* [Following Lagrange's role we may write the second term of jS in the form 

The quantity outside the integral sign is to be taken between the given limits 
and is zero, since vanishes at each limit. Let us now put 

It is clear that this value of X cannot be infinite between the limits of integra- 
tion unless u vanishes. For by hypothesis u and u' are both finite and the co- 
efficients £ii and C^ in the Lagrangian function are also finite. It then eai^y 
follows from the equation 

that (Cji + 2X)» = Bii {A^ + 2X'). 

Hence the second term of S becomes 



lJBn{^-eiy.u, 



which is the result in the text. This might also have been deduced from Jacobi's 
general transformation with one independent variable given in Prof. Jellett's 
Calculus of Variations or Prot Price's Differential Calculus, 

li u=^0 between the limits of integration this transformation fails, but, as is 
evident from the argument in the text, we choose ^=m to represent a neighbouring 
free motion such that u=0 just before the system reaches A ; also B is so placed 
that the next instant at which u=0 is after the system has passed B.] 
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14. [We shall conclude the chapter with the application of 
this criterion of stability to some simple case. 

A particle describes a circular orbit about a centre of force 
situated in the centre. It is required to deduce the conditions of 
stability as to space from the Characteristic fum>tion V. 

Let a be the radius of the circle, n the angular velocity of the 
particle about the centre 0. Let ^(a)be the law of force. Let 
A and B be two points taken on the circular orbit, and let the 
particle be conducted from -4 to 5 by some neighbouring path 
with the same energy as in the circular orbit. Let r=a-\- phe 
th6 radius vector of this path, corresponding to any angle 0. 

If V be the velocity at any point of this path, we easily find 



«=«»{' -^(1^^')^- 



If « be the arc of the path, we have 

ds ^ ^ 1 /dp\' 

If the angle A OB = ^, we therefore have 

F = jvds = aV/3 + |£|(g) -/p j d0, 

where 1)' = '^ + 3. 

If the neighbouring path be a free path described with the 
same energy, its equation is 

p = 2/ sin pd, 

where is measured from the radius vector OA, and L is an 
arbitrary constant. This free path will cut the circle again in 
some point C. If the angle -4 0(7= 7, we have j>7 = tt. . 

If B coincide with C, we find by substituting this value of /o in 
the expression for V, that the second term of Fis zero. If -B be 
beyond G but such that the angle COB is less than 7, draw two 
free paths one from A as before and the other backwards from B 
to meet the former in some point P. Then the angle AOP^\p. 
If the particle be conducted from A to P along one path and 
from F to B along the other, we find that the excess of the 
action over the action in the circular arc AB is equal to 

—-L^sinp^, 

Since pyS is greater than tt and less than 27r, this excess is 
negative. The action along the circular arc is therefore not a 
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minimum if jS be beyond the first intersection of a neighbouring 
free path. 

Lastly we may show that the action is a minimum if B lie 
between A and C. To prove this we write the integral in the 
expression for V in the form 

The first of these two terms is zero since p vanishes at each 
limit Following Lagrange's rule we make 

The integral then becomes 

This is always positive and the action along the circular arc is 
a minimum. The argument however requires that \ should not 
become infinite between the limits of integration. It is easy to 
see that 

where jE is a constant to be chosen at our pleasure. But if fi 
exceed - it is impossible to choose JE^ so as to keep \ finite be- 
tween the limits of integration. Hence the action is a minimum 
only if the angle A OB subtended by the limiting positions at the 

centre is less than -. The circular orbit is therefore stable if p* 

. , P 

be positive. 

If, however, p^ be negative, the expression for \ changes its 
character. Writing — g' for p^, we find 

thus the value of \ can be chosen so as not to be infinite for all 
positive values of 0. In this case the function F is a minimum 
for all arcs AB however distant B may be from A. The circular 
motion is therefore unstable.] 
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